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ABSTRACT 

We present new Keck/DEIMOS spectroscopic observations of hundreds of individual stars along 
the sightline to Andromeda's first three discovered dwarf spheroidal galaxies (dSphs) - And I, II, and 
III, and leverage recent observations by our team of three additional dSphs, And VII, X, and XIV, 
as a part of the SPLASH Survey (Spectroscopic and Photometric Landscape of Andromeda's Stellar 
Halo). Member stars of each dSph are isolated from foreground Milky Way dwarf and M31 field 
contamination using a variety of photometric and spectroscopic diagnostics. Our final spectroscopic 
sample of member stars in each dSph, for which we measure accurate radial velocities with a median 
uncertainty (random plus systematic errors) of 4 - 5 km s _1 , includes 80 red giants in And I, 95 in 
And II, 43 in And III, 18 in And VII, 22 in And X, and 38 in And XIV. The sample of confirmed 
members in the six dSphs are used to derive each system's mean radial velocity, intrinsic central 
velocity dispersion, mean abundance, abundance spread, and dynamical mass. 

This combined data set presents us with a unique opportunity to perform the first systematic 
comparison of the global properties (e.g., metallicities, sizes, and dark matter masses) of one-third of 
Andromeda's total known dSph population with Milky Way counterparts of the same luminosity. Our 
overall comparisons indicate that the family of dSphs in these two hosts have both similarities and 
differences. For example, we find that the luminosity - metallicity relation is very similar between L ~ 
10 5 - 10 7 Lq, suggesting that the chemical evolution histories of each group of dSphs is similar. The 
lowest luminosity M31 dSphs appear to deviate from the relation, possibly suggesting tidal stripping. 
Previous observations have noted that the sizes of M31's brightest dSphs are systematically larger 
than Milky Way satellites of similar luminosity. At lower luminosities between L = 10 4 - 10 6 Lq, 
we find that the sizes of dSphs in the two hosts significantly overlap and that four of the faintest 
M31 dSphs are smaller than Milky Way counterparts. The first dynamical mass measurements of six 
M31 dSphs over a large range in luminosity indicates similar mass-to-light ratios compared to Milky 
Way dSphs among the brighter satellites, and smaller mass-to-light ratios among the fainter satellites. 
Combined with their similar or larger sizes at these luminosities, these results hint that the M31 
dSphs are systematically less dense than Milky Way dSphs. The implications of these similarities and 
differences for general understanding of galaxy formation and evolution is summarized. 

Subject headings: dark matter — galaxies: dwarf — galaxies: abundances — galaxies: individual 
(And I, And II, And III, And VII, And X, And XIV) — galaxies: kinematics - 
galaxies: structure — techniques: spectroscopic 
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1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies represent important 
laboratories for constraining feedback processes in galax- 
ies and for testing cosmological models on the smallest 
scales (Klypin et al. 1999; Moore et al. 1999; Bullock, 
Kravtsov, & Weinberg 2000). These low luminosity sys- 
tems (—1.5 < My < —13) show no evidence of ongo- 
ing star formation, contain little or no interstellar mat- 
ter, and are almost always distributed around massive 
hosts, such as the Milky Way. Within the framework of 
the ACDM paradigm, theories suggest that the observed 
dSph satellites are the present-day counterparts to pre- 
viously accreted objects that were tidally destroyed to 
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build the spheroids and stellar halos of massive galaxies 
(Bullock et al. 2001; Bullock & Johnston 2005). This 
process has been directly verified in the outskirts of gi- 
ant galaxies with the discovery of luminous substructure 
in the form of tidal streams. For example, the accre- 
tion mechanism, where these dwarf galaxies are shredded 
and eventually melded into diffuse stellar components, is 
abundant in the Milky Way (e.g., the Sagittarius, Or- 
phan, and Cetus Polar Streams - Ibata, Gilmore, & Irwin 
1994; Belokurov et al. 2007; Newberg, Yanny, & Willet 
2009), in the Andromeda (M31) spiral galaxy (e.g., the 
Giant Southern Stream - Ibata et al. 2001), as well as 
in several other systems such as NGC 5907 (Martinez- 
Delgado et al. 2008) and NGC 4013 (Martinez-Delgado 
et al. 2009). 

Despite representing the most common morphological 
class among all galaxies in the Universe, the overall count 
of known dSph satellites around the Milky Way and M31 
are significantly lower than first-order expectations. Sim- 
ulations predict large galaxy halos to host thousands of 
dark matter substructures with masses similar to those 
of dSph galaxies (Klypin et al. 1999; Diemand, Kuhlen, 
& Madau 2007; Springel et al. 2008). While models ex- 
ist for explaining the mismatch, they can only be tested 
with precise mass measurements of the observed dwarfs. 
Recently, large observational efforts have successfully un- 
covered many new dSphs around both the Milky Way 
(e.g., Willman et al. 2005; Zucker et al. 2006a; 2006b; 
Sakamoto & Hasegawa 2006; Belokurov et al. 2006; Be- 
lokurov et al. 2007; Irwin et al. 2007; Walsh, Jerjen, & 
Willman 2007; Liu et al. 2008; Belokurov et al. 2008; Be- 
lokurov et al. 2009) and M31 (Martin et al. 2006; Zucker 
et al. 2007; Ibata et al. 2007; Majewski et al. 2007; Irwin 
et al. 2008; McConnachie et al. 2008; Martin et al. 2009), 
effectively doubling the total sample size in the past few 
years. Follow up spectroscopic observations of individ- 
ual stars in the new Milky Way satellites (e.g., by Kleyna 
et al. 2005; Munoz et al. 2006; Simon & Geha 2007, Mar- 
tin et al. 2007; Koch et al. 2009; Geha et al. 2009) have 
established that these systems are the most dark matter 
dominated objects known, and also provided a means to 
test proposed solutions of the "missing satellites prob- 
lem" (see also Strigari et al. 2008a; 2008b). Unfortu- 
nately, it is difficult to draw general conclusions based 
on the Milky Way system alone. 

Recently, a number of large surveys have targeted the 
outskirts of M31 with wide field imagers on 4-meter class 
telescopes (e.g., Ferguson et al. 2002; Osthcimcr 2003; 
Ibata et al. 2007; McConnachie et al. 2009). These stud- 
ies, combined with spectroscopic follow up such as the 
SPLASH Survey (Spectroscopic and Photometric Land- 
scape of Andromeda's Stellar Halo), have resulted in 
ground-breaking discoveries related to the properties of 
M31's stellar halo, which itself was discovered in the 
surveys (Guhathakurta et al. 2005; Irwin et al. 2005). 
For example, the field halo of M31 exhibits a power- 
law surface brightness profile similar to the Milky Way 
(Osthcimer 2003; Guhathakurta et al. 2005; Irwin et al. 
2005), is filled with substructure (Ibata et al. 2007; 
Gilbert et al. 2007, 2009), and is metal-poor in its out- 
skirts (Kalirai et al. 2006a; Chapman et al. 2006; Koch 
et al. 2008). These recent surveys have also led to the dis- 
covery of new dSph galaxies in M31, over 35 years after 
van den Bergh (1972) found the first such systems. Sim- 



ilar to the discovery rate in the Milky Way, the past five 
years have led to thirteen new galaxies, And IX (Zucker 
et al. 2004), And X (Zucker et al. 2007), And XI-XIII 
(Martin et al. 2006), And XIV (Majewski et al. 2007), 
And XV-XVI (Ibata et al. 2007), And XVII (Irwin et al. 
2008), And XVIII-XX (McConnachie et al. 2008), and 
And XXI-XXII (Martin et al. 2009). 11 The total pop- 
ulation of known dSphs is therefore similar in M31 and 
the Milky Way, and considering the spatial coverage of 
the surveys in both galaxies (e.g., the PAandAS Survey - 
Ibata et al. 2007; McConnachie et al. 2008; McConnachie 
et al. 2009), the overall census is incomplete by at least 
a factor of a few. 

Given the distance to M31, 780 kpc, most of our infor- 
mation on the properties of its dSph population has come 
about from studies of the bright red giant branch (RGB) 
stars in these systems. For example, McConnachie & 
Irwin (2006) present a detailed analysis of the struc- 
tural properties of six of these galaxies and find both 
similarities and differences with the Galactic population. 
For example, similar to the Galactic system, M31 dSphs 
with higher central surface brightnesses are found fur- 
ther from their host. A key difference however is that 
the scale radii of the M31 dSphs is approximately twice 
as large as the Milky Way dSphs at the same luminosity. 
McConnachie & Irwin (2006) also show that the central 
surface brightness of the M31 satellites is smaller than 
that of comparable Milky Way systems. The only Hubble 
Space Telescope HST observations of the M31 satellites 
are the WFPC2 studies of Da Costa et al. (1996; 2000; 
2002), who obtained optical photometry of stars in each 
of And I, II, and III down to the level of the horizontal 
branch. Their study provides the first constraints on the 
age and abundance distributions of stars in these three 
systems. Indeed, Da Costa et al. find that, similar to 
the Milky Way dSphs, the M31 satellites have also had 
diverse evolutionary histories. 

Unlike for the Milky Way satellites, there exist very 
few spectroscopic studies of the M31 dSph system, most 
of which were only able to characterize a few stars in each 
galaxy. Cote et al. (1999) obtained the first such spec- 
tra using HIRES on Keck. They measured velocities for 
seven RGB stars in And II and estimated both the radial 
velocity of the galaxy and the central velocity dispersion 
from this small sample. Guhathakurta, Reitzel, & Grebel 
(2000) also obtained spectra of a few tens of individual 
RGB stars in each of And I, III, V, and VI. Their study 
led to a dynamical mass estimate of M31 based on the 
mean radial velocity of each dSph (Evans et al. 2000). 
Unfortunately, the data set could not be used to probe 
the internal kinematics of the satellites given low num- 
bers of stars in each galaxy, and large uncertainties in 
the velocity measurements. Guhathakurta, Reitzel, & 
Grebel (2000) also looked at several stars in And VII 
with HIRES on Keck and measured both the radial ve- 
locity and the total velocity dispersion of the galaxy (e.g., 
the combined true intrinsic dispersion plus the dispersion 
due to velocity errors). Recently, Chapman et al. (2005) 
and Chapman et al. (2007) have used Keck/DEIMOS to 
spectroscopically confirm six to eight members of And IX 
and And XII, however these data do not constrain the 
internal kinematics of the dSph. For example, the inclu- 

11 And XVIII is a distant background Local Group dSph. 
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TABLE 1 



Date 


Mask Pointing center: 


Field PA 


Exp. 


No. Sci. 




",72000 5.12000 


(°E of N) 


Time 


Targets 1 




^h.m.s^ ^o.f.ff^ 




(s) 




2005 Nov 05 dl-1 (And I) 00:45:48.61 +38:05:46.8 


+0.0 


3 x 1,200 


153 


2006 Sep 16 


dl_2 (And I) 00:46:13.95 +38:00:27.0 


+90.0 


3 x 1,200 


150 


2005 Sep 06 


d2_l (And II) 01:17:07.46 +33:29:25.1 


-90.0 


3 x 1,200 


139 


2005 Sep 06 


d2_2 (And II) 01:16:43.29 +33:34:25.8 


+0.0 


3 x 1,200 


139 


2005 Sep 08 


d3.1 (And III) 00:36:03.83 +36:27:27.4 


+90.0 


3 x 1,200 


128 


2005 Sep 08 


d3.2 (And III) 00:35:39.61 +36:21:41.8 


+0.0 


3 x 1,200 


126 



1 10 of the stars in And I, 17 in And II, and 15 in And III were 
observed on both masks. 



sion of one marginal member in the outer radial bin of 
And IX inflates the velocity dispersion by a factor of two. 
A re-analysis of these data by Collins et al. (2009) estab- 
lishes upper limits on the velocity dispersions of both 
galaxies 12 , and they also do not resolve the dispersions 
of And XI and XIII. Finally, Letarte et al. (2009) have 
presented a new analysis of seven and eight spectroscop- 
ically confirmed stars in And XV and And XVI. The 
radial velocity uncertainties in these measurements are 
large, ranging from a v = 6 - 25 km s _1 , with one star 
as high as a v = 52 km s _1 . Summarizing, of the entire 
18 galaxies comprising the known M31 dSph population, 
And II and And VII are the only systems with secure ve- 
locity dispersion measurements, based on high resolution 
studies of 7 and ^20 stars each. 

In this paper, we present the first results from a large 
project aimed at spectroscopically observing hundreds 
of stars in the M31 dSph system, as a part of the over- 
all SPLASH Survey. These data will allow the first de- 
tailed comparison of the properties of M31 satellites to 
the Milky Way dSph system. For example, we will ac- 
curately measure velocities of (a minimum of) dozens of 
individual stars in each galaxy, thereby resolving their 
intrinsic central velocity dispersions and leading to esti- 
mates of the total system masses. For the best data sets 
(>500 velocities in a single dSph), we will also probe for 
any radial changes in the velocity dispersion, and test 
for rotation, thereby constraining dynamical models for 
dSphs (e.g., Kroupa 1997). Such high quality data sets 
can also be used to better understand how luminosity 
maps to dark matter subhalo mass. 

In addition to the kinematical analysis, spectroscopic 
observations of large numbers of stars in each dSph can 
provide a detailed analysis of radial trends in the chemi- 
cal abundances of stars, and perhaps even constrain the 
a-abundance of these satellites. Therefore, these ob- 
servations directly relate to a better understanding of 
both population gradients and the evolutionary history 
of dSphs. Comparisons of the results from the survey 
with the well studied dSphs of the Milky Way will aid in 
our understanding of whether the observed differences in 
the satellite systems (see § 6) are caused by in situ pro- 
cesses within the galaxies, or whether they are shaped 
by different interaction histories of the dSphs with their 

12 At the time of the writing of this paper, the Collins et al. 
(2009) study appears as a preprint. The l-cr error bars on the 
velocity dispersions are consistent with c v = km s for both 
dSphs. 



massive hosts. 

In this first paper, we address the internal kinematics, 
chemical compositions, sizes, and total masses of And I, 
II, III, VII, X, and XIV, based on a large sample size of 
individual radial velocities. Future papers will address 
other M31 dSphs, for which we have already collected a 
substantial amount of Keck/DEIMOS observations. The 
layout of this paper is as follows. We describe the pho- 
tometric and spectroscopic observations of these dSphs 
as well as the data reduction in the next section. This 
includes a detailed discussion of the uncertainties in the 
velocity measurements, which are critical to understand 
in order to resolve the intrinsic dispersions of the (kinc- 
matically cold) satellites. Membership of stars in the 
satellites is established and velocity histograms of the 
confirmed red giants are presented in § 3. The kinematics 
are analyzed to yield mass-to-light ratios for the galax- 
ies in §4. Next, we measure the chemical abundance of 
each dSph in § 5. Unlike previous studies, we measure 
the abundances of each confirmed dSph star using two 
independent methods, photometrically and spectroscop- 
ically. The results from this overall study are discussed 
in § 6, where we compare the properties of all six M31 
dSphs (e.g., their sizes, metallicities, and masses) with 
measured trends for Milky Way satellites. The main re- 
sults from the paper are summarized in § 7. 

2. THE DATA SET 

The imaging and spectroscopic observations for three 
of the six dSphs in this study are being presented for 
the first time (And I, II, and III), whereas the data sets 
for And VII, X, and XIV have been (at least partially) 
presented in previous studies. For these three galaxies, 
however, we present a new analysis of the spectra in the 
present paper. We begin by first introducing the new 
data for And I, II, and III, followed by a short summary 
of the previously presented data for And VII, X, and 
XIV, including any relevant new analysis. 

2.1. Photometric and Spectroscopic Observations - 
And I, II, and III 

We imaged each of And I, II, and III with the wide-field 
Mosaic camera on the Kitt Peak National Observatory 
(KPNO) 4-m telescope from 1998 to 2002 (Ostheimer 
2003). This camera subtends an angular size of 36' x 36', 
much larger than the King limiting diameter of And I (r/ 
= 10.4') and And III (r, = 7.2'), and slightly smaller 
than that of And II (r; = 22.0'). The observations were 
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Fig. 1. — The distribution of stars from the KPNO wide-field Mosaic observations of each dSph are shown as small points. The selection 
function used to isolate these stars involved a simple cut on the color-magnitude diagram (CMD) to enhance the dSph RGB with respect 
to the field population, which is still abundantly present in the maps. Superimposed on the photometric data are the objects targeted with 
Keck/DEIMOS spectroscopy (darker points). 



obtained in the Washington System (M, T 2 ) bands and 
converted to Johnson-Cousins (V, /) using the relations 
given in Majewski et al. (2000). The exposure times were 
set at 900 seconds in the M filter and 3,600 seconds in 
the T filter, reaching a photometric depth of I ~ 24. 
The most luminous (metal-poor) RGB stars have intrin- 
sic magnitudes of Mj ~ —4, and therefore at the distance 
of M31 (780 kpc, (m-M)o = 24.5), these stars will have 
la = 20.5. Our imaging observations are therefore well 
suited to probe several magnitudes of the RGB in each 
of the dSphs. 

We also obtained imaging observations of the same 
fields using the intermediate- width DDOdl filter (Ma- 
jewski et al. 2000). This filter is centered on the surface- 
gravity sensitive Mgb and MgH stellar absorption fea- 
tures, and therefore provides an efficient means to dis- 
criminate foreground Milky Way contaminants from M31 
dSph RGB stars (Ostheimer 2003; Gilbert et al. 2006; 
Guhathakurta et al. 2006). The filter is especially effi- 
cient at boosting our odds of targeting true dSph stars 
that are located far from the core of each galaxy, where 
the surface brightness of members if low relative to 
the foreground veil of Milky Way dwarf stars. These 
stars, in the outskirts of the dSphs, represent key data 
points when studying radial trends. The combination of 
the wide-field broadband imaging observations with co- 
spatial DDOhl filter exposures provides the necessary 
framework to conduct the first largely uncontaminated 
study of the properties of M31 dSphs. 

With this unique photometric data set in hand, we 
spectroscopically targeted the brightest RGB candidates 
in each of And I, II, and III with the DEIMOS multi- 
object spectrograph (Faber et al. 2003) on the Keck II 
10-meter telescope. To do this, we first measured the po- 
sitions of all stars in the KPNO Mosaic images and trans- 
formed the positions into equatorial coordinates using 
the USNO-B Guide Star Catalogue as a reference (Monet 
et al. 2003). These positions were then mapped onto the 
DEIMOS mask. The spectroscopic data were collected 
in late 2005, as summarized in Table 1. The observations 
were obtained using the standard setup described in de- 
tail in several previous papers (e.g., Gilbert et al. 2006, 



Guhathakurta et al. 2006 and Kalirai et al. 2006b). To 
summarize a few key details, the spectroscopic observa- 
tions were obtained using the 1200 lines mm -1 grating 
at a central wavelength of 7800 A. The resolution of the 
spectra for a typical seeing of 0'.'8 FWHM is 1.3 A and 
the spectral range is «6400 - 9100 A, depending exactly 
on the position of each object relative to the mask center. 

In total, we measured a spectrum for 835 objects in 
six masks (two per galaxy), where each mask subtends 
roughly 16' x 4'. Of these, 42 objects were measured 
on multiple masks to assess our radial velocity accuracy 
(see §2.3). The targets were chosen over a magnitude 
range extending from just above the tip of the RGB down 
to ^2 magnitudes below the tip in each satellite. The 
mean spectral S/N in the entire sample of observations 
is ^5 per pixel. The DD051 filter observations were 
used to screen objects in the selection process to ensure 
that high probability RGB stars were targeted before 
lower probability objects (which turn out to be mostly 
foreground Milky Way dwarfs). 

The positions and orientations of the Keck/DEIMOS 
spectroscopic fields relative to And I, II, and III were 
carefully placed to allow both a study of stars near the 
center of each dSph and also to probe radially outward 
from the core. To illustrate these pointings, we present 
starcount maps of each galaxy taken from our wide-field 
KPNO Mosaic imaging data in Figure 1. For each galaxy, 
we make a very rough cut on the photometric data to iso- 
late stars located along the RGB in the color-magnitude 
diagrams (CMDs), which are presented later in §3. The 
spectroscopic targets are shown as darker points. 

2.1.1. Existing Observations of And VII, X, and XIV 

The data sets for And VII, X, and XIV differ from 
the three dSphs discussed above in several ways, and 
so we briefly summarize this sample here. First, for 
And VII, the photometry was obtained by Grebel & 
Guhathakurta (1999) using Keck/LRIS. The galaxy is 
one of the most luminous M31 satellites known, with an 
integrated brightness of My = —13.3 ± 0.3, and so sev- 
eral of the brightest red giants in the photometric sample 
were targeted with Keck/HIRES for spectroscopy. These 
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observations utilized special purpose multislit masks with 
a half dozen slitlets each, with each slitlet having a length 
of 1 - 2'.' 5. The technique employed is described in de- 
tail in Sneden et al. (2004), in the context of similar 
observations of the globular cluster M3. Note, Grebel 
& Guhathakurta (1999) refer to this galaxy as the Cas- 
siopeia dSph in their paper. 

The second galaxy in this existing sample is And X, 
a newly discovered low-luminosity M31 dSph from the 
Sloan Digital Sky Survey (Zucker et al. 2007). We used 
the imaging data set described by Zucker et al. (2007) to 
target two Keck/DEIMOS spectroscopic masks in this 
galaxy, using the same setup as discussed above (and 
below) for And I, II, III. The details of the observations, 
data reduction, and results are presented in Kalirai et al. 
(2009). 

The final galaxy in the sample is And XIV, a system 
that we discovered with KPNO/Mosaic three years ago. 
This galaxy is one of the most remote systems known 
in M31, located at a projected distance of 162 kpc and 
a line of sight distance that places it behind M31 by 
~80 kpc. We targeted And XIV with two Keck/DEIMOS 
spectroscopic masks, again, using an identical setup as 
discussed above (and below) for And I, II, and III. The 
details of the observations, data reduction, and results 
for this dSph are presented in (Majewski et al. 2007). 

2.2. Data Reduction 

Data reduction of the six DEIMOS multi-slit masks 
being presented for the first time in the present study 
(And I, II, and III) are based on the spec2d software 
pipeline (version 1.1.4) developed by the DEEP2 Galaxy 
Rcdshift Survey team at the University of California- 
Berkeley for that project (see Faber et al. 2007 for more 
information). Briefly, internal quartz flat- field exposures 
are obtained and used to rectify the curved raw spectra 
into rectangular arrays by applying small shifts and in- 
terpolating in the spatial direction. A one-dimensional 
slit function correction and two-dimensional flat-field and 
fringing correction are applied to each slitlet. Using 
the DEIMOS optical model as a starting point, a two- 
dimensional wavelength solution is determined from mul- 
tiple Kr-Ar-Ne-Xe arc lamp exposures with residuals of 
order 0.01 A. Each slitlet is then sky-subtracted exposure 
by exposure using a B-splinc model for the sky. The indi- 
vidual exposures of the slitlet are averaged with cosmic- 
ray rejection and inverse- variance weighting. Finally one 
dimensional spectra are extracted for all science targets 
using the optimal scheme of Home (1986) and are re- 
binned into logarithmic wavelength bins with 15 km s _1 
per pixel. We note that the modifications to this pipeline 
discussed by Simon & Geha (2007) and Gilbert et al. 
(2009) were used in the present reductions. Specifically, 
the cosmic ray rejection algorithm was altered to allow 
alignment of individual two dimensional exposures in the 
spatial direction before co-addition. 

In Figure 2, we present several of our final 1- 
dimensional spectra for stars in each of the three M31 
dSphs. Only a 500 A section of the spectrum centered 
near the Can triplet is shown in each case for the sake 
of clarity. These absorption lines, as well as other spec- 
tral features, are clearly visible in the data for bright and 
faint objects in each galaxy. 



2.3. Radial Velocity Measurements 

Radial velocities are measured for all extracted one- 
dimensional spectra by cross-correlating the observa- 
tions with a series of high signal-to-noise stellar tem- 
plates. The stellar templates were observed using the 
same DEIMOS setup described above, and range in spec- 
tral type from F8 III to M8 III, and also include subgiants 
and dwarfs. For the faintest objects with lower signal- 
to-noise, this cross-correlation does not yield meaningful 
results and therefore we manually check all of the tem- 
plate fits and only include those spectra in which the 
template matched at least two different spectral features 
(see Guhathakurta et al. 2006 for more details). Of the 
initial sample of 835 spectra, this left 426 spectra after 
eliminating galaxies, 148 in And I, 136 in And II, and 
142 in And III. These numbers include 1 duplicate in 
the And I observations, 1 in And II, and 8 in And III. 
The duplicates represent observations of the same star 
on both masks, where both spectra yielded a reliable ve- 
locity measurement from the cross correlation. 

Below the set of sample spectra for each dSph in Fig- 
ure 2 is shown the inverse of the variance spectrum (red), 
computed under the idealized assumption that it is dom- 
inated by Poisson statistics in the sky+object counts. 
The spectral measurements presented in this paper — the 
cross-correlation analysis used to measure radial veloc- 
ities and the measurement of absorption line strengths 
(see § 5.2) — properly take into account this dependence 
of the uncertainty in the spectral flux on wavelength. 
The dips in the inverse variance spectrum correspond to 
bright night sky emission lines, shifted in wavelength to 
account for the shift in the science spectra from the ob- 
server frame to the rest frame. Occasionally, systematic 
errors in sky subtraction can cause the actual uncertainty 
to be greater than the Poisson errors shown here. The 
inverse variance spectra are plotted to help the reader 
distinguish between real features (e.g., the Can triplet 
absorption lines) and noise artifacts for the fainter stars 
in our sample. 

The uncertainties in our velocity measurements are 
empirically estimated to be ^10 km s , based on the 
comparison of independent observations from these du- 
plicates. While this uncertainty is many times less than 
the velocity dispersion in bulge and halo- like populations, 
it is comparable to, if not larger than, the expected ve- 
locity dispersion in small dSphs. A significant fraction 
of the radial velocity uncertainty budget derives from 
random and systematic astrometric errors. As discussed 
earlier, in designing our spectroscopic masks we used the 
USNO-B Guide Star Catalogue to calculate absolute po- 
sitions of all stars. The internal accuracy of this cata- 
logue is known to be ~0"2 (Monet et al. 2003). At the 
DEIMOS pixel scale (0'.'1185 per pixel), this uncertainty 
alone translates to 1.7 pixels. Civen that the disper- 
sion of the 1200 lines mm -1 grating is 0.33 A pixel -1 
and the anamorphic factor is 0.6 at 8000 A, a 1.7 pixel 
uncertainty results in a 0.3 A wavelength uncertainty. 
Therefore, a random misalignment of the centroid of a 
star by 0'.'2 in the slit can lead to a velocity error of 
12 km s _1 . The positional accuracy of a given target 
will also depend on the overall solution of the astromet- 
ric transformation and therefore depends on the density 
of astrometric reference stars available at a particular 
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Fig. 2. — Five sample spectra of RGB stars in And I, II, and III corrected to zero velocity. To illustrate the quality of the full data set, 
these spectra have been chosen over a magnitude range extending from the brightest to the faintest stars in our spectroscopic sample. The 
wavelength window of 500 A has been centered near the Can triplet absorption lines (8498, 8542, and 8662 A) which are clearly detected 
in all objects. The red curve curve shows a sample inverse variance spectrum for each dSph, calculated under the assumption that the 
errors are dominated by Poisson statistics in the sky+object counts, and intended to help distinguish between real features (e.g., Can 
triplet absorption lines) and noise artifacts for the faintest stars plotted. The full spectral coverage of these data is ~2,700 A, which allows 
multiple features to be used in the cross-correlation to yield accurate radial velocities. The spectra have been smoothed using a 7 pixel 
boxcar function for illustration purposes only. 



location in the imaging frame. Finally, we note that a 
number of other effects can also lead to uncertainties in 
the astrometry. These can include positional uncertain- 
ties in the images (e.g., if one or more reference stars 
have an appreciable proper motion), slight imperfections 
in the mechanical cutting process that produces the slits, 
and even small misalignments (i.e., rotations or shifts) of 
the entire mask during the observations. 

Fortunately, these astrometric uncertainties leave a sig- 
nature that can be used to calculate the amount that the 
radial velocity measurement is offset from the true value. 
This method, as discussed by Sohn et al. (2007) and Si- 
mon & Geha (2007), utilizes atmospheric absorption fea- 
tures that are superimposed on the stellar spectra. The 
strongest of these telluric lines are in the A-band (7600 
- 7630 A). To correct for the mis-centering, we observed 
several bright, rapidly rotating, hot standard stars with 
few emission or absorption lines (e.g., the Wolf Rayet 
star Wolf 1346). The near featureless spectra of these 



stars displays a very clean detection of the A-band ab- 
sorption. The telluric templates were created by allowing 
these stars to drift perpendicularly across the slit during 
each short exposure, ensuring that the starlight evenly 
fills the slit. Next, we averaged the spectra for all such 
standards and cross-correlated a windowed region of the 
coadded spectra at the wavelength of the A-band absorp- 
tion to the observed stellar spectra. For a star that is 
perfectly centered in the slit, this cross-correlation pro- 
duces a near perfect match indicating that no offset is 
required. However, for a star that is mis-centered in the 
slit, the observed wavelength of the A-band is offset rel- 
ative to the standard star template. We use this offset 
to calculate a velocity shift and apply this correction to 
each of our 426 radial velocities. The mean offset for 
each mask ranged from 0.7 - 10.2 km s _1 , with a stan- 
dard deviation of 5 - 13 km s _1 . The final addition to 
each of our velocity measurements includes a heliocen- 
tric correction to convert our geocentric velocities into a 
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Fig. 3. — Repeat measurements of the same stars in our And I, II, 
III, and X data sets are used to define a normalized error, <tjv equal 
to the ratio of the difference in velocity of the duplicates and the 
quadrature sum of all error terms. In addition to the Monte Carlo 
error estimates, we find that an additional error (e) is needed to 
reproduce a unit Gaussian in the resulting distribution of crjv (solid 
histogram). As discussed in § 2.4, we formally add e = 2.2 km s 
in quadrature to each of our radial velocity errors to match the 
results from Simon & Geha (2007), which are based on a much 
larger sample of repeat measurements taken with a similar setup 
to this work. 



heliocentric frame of reference. 

2.4. Radial Velocity Uncertainties 

The observed velocity dispersions of the M31 satellites 
in this study reflect a combination of the actual intrinsic 
dispersions, plus the velocity errors. The known intrin- 
sic dispersion of other Local Group satellites are often 
5-10 km s _1 , which is comparable to the expected er- 
rors in our velocity measurements. It is therefore very 
important to have a good understanding of the true un- 
certainties in our measurements. We measure these un- 
certainties using a Monte Carlo method, where we first 
take each stellar spectrum and add noise to each pixel 
scaled by the estimated variance in that pixel (also see 
discussion in Simon & Geha 2007). This procedure is 
repeated 1000 times assuming the variance in each pixel 
is distributed according to Poisson statistics, and the ve- 
locity and telluric correction is re-calculated after each 
run of the simulation. The error in the velocity (cr) is de- 
fined as the square root of the variance in the recovered 
mean velocity. 

We next compare this error measurement to the known 
uncertainty, measured as the difference in velocity be- 
tween independent measurements of the same star (e.g., 
v\ and V2). For this, we use the 10 duplicates with re- 
liable velocity measurements in this data set, as well as 
the 7 duplicates in our And X data set (Kalirai et al. 
2009). For these stars, we define a normalized error, cat, 
as the ratio of the velocity difference between duplicate 
measurements and the quadrature sum of all error contri- 



FlG. 4. — The total uncertainty in our velocity measurements 
for all stars along the sightline to And I, II, and III is shown as 
a function of the mean spectral S/N per pixel. The uncertainties 
reflect both random and systematic velocity error contributions, 
the two being added in quadrature as discussed in § 2.4. Even for 
spectra with S/N = 5, we are able to measure velocities accurately 
to a v = 4 km s —1 . 



butions (<7i , 1T2 , and \/2e) . The latter includes the Monte 
Carlo errors in each of the two velocity measurements, as 
well as an additional term representing any other error 
that we have not accounted for. The distribution of o~n 
for all duplicates is shown in Figure 3, and in order to re- 
produce a unit Gaussian, we measure e to be <3 km s . 
For our final velocity errors, we choose to formally add an 
e of 2.2 km s _1 in quadrature with the Monte Carlo er- 
rors. Although our measured value is slightly larger than 
this, 2.2 km s _1 reflects the e measured by Simon & Geha 
(2007) from a similar Keck/DEIMOS study, but with a 
factor of three larger data set of duplicate measurements. 
With only 17 data points, our measured cfn is much more 
susceptible to a few outliers than their sample. We note 
that the resulting velocity errors from this analysis, and 
with this value of e, have been directly tested by Simon 
& Geha (2007) by comparing to high dispersion observa- 
tions of a Milky Way globular cluster and dSph galaxy, 
each with velocity uncertainties of <1 km s _1 . The re- 
sults from the Keck/DEIMOS analysis are in excellent 
agreement (<lc) with the high resolution studies. 

In Figure 4 we present the final error distribution of all 
velocity measurements in our data set, as a function of 
the spectral S/N (per pixel). As expected, the highest 
S/N stars have a velocity error equal to the "floor" es- 
tablished above (2.2 km s _1 ) whereas the error increases 
to ~4 km s _1 for stars with S/N — 5, and quickly in- 
creases for lower S/N spectra. The median uncertainty 
across the full sample is 4.4 km s _1 . 

We stress that the analysis presented above in carefully 
understanding the velocity errors in these data is essen- 
tial to accurately probe the internal kinematics of dwarf 
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Fig. 5. — The velocity histograms of And VII (left - Guhathakurta, Reitzel, & Grebel 2000), And X (middle - Kalirai et al. 2009), and 
And XIV (right - Majcwski et al. 2007). The properties of each of these satellites, as measured from these kinematical data, are discussed 
in the text. 



galaxies. The uncertainty introduced by even very small 
astrometric errors (e.g., at the 0'.'05 - O'.'l level) in the 
centroids of targets leads to radial velocity uncertainties 
that are larger than the dispersions of some satellites. 
Furthermore, these uncertainties, for a given spectrum, 
systematically offset the wavelengths of all spectral fea- 
tures by the same amount. The rms scatter between 
velocity measurements, as measured from different ab- 
sorption lines, can therefore not provide an accurate as- 
sessment of the overall uncertainty in radial velocity. 

2.5. Radial Velocity Measurements for 
And VII, X, and XIV 

The method used to reduce the And X and XIV 
Keck/DEIMOS spectroscopic data, including the mea- 
surement of radial velocity and uncertainty, is very sim- 
ilar to the description above and has been presented 
in Kalirai et al. 2009 and Majewski et al. 2007). For 
And VII, the Kcck/HIRES data were analyzed slightly 
differently. The sky subtraction was performed by com- 
bining the few pixels of sky in all of the slits into a master 
sky spectrum. The velocities were extracted by cross cor- 
relating the And VII target spectra with radial velocity 
standard stars. The latter were observed through the 
same slitlets on the same HIRES masks as the science 
observations. Although the S/N of the spectra are quite 
low, the wide baseline in wavelength allowed reliable ve- 
locities to be measured for 18 of the 21 And VII RGB 
candidates. The uncertainties in the individual velocity 
measurements are <1.5 km s . 

3. ESTABLISHING MEMBER STARS IN EACH DSPH 

Membership of stars in the three previously observed 
dSphs, And VII, X, and XIV, have been discussed in 
Guhathakurta, Reitzel, & Grebel (2000), Kalirai et al. 
(2009), and Majewski et al. (2000), and so we refer to 
those papers for the details (the analysis of And X and 
XIV is identical to the discussion below). Summarizing, 
accurate radial velocities and uncertainties were estab- 
lished for 18 member stars in And VII, 22 member stars 
in And X, and 38 member stars in And XIV. These pa- 



pers, as well as Grebel & Guhathakurta (1999), present 
the CMDs for each of these satellites. We summarize the 
three radial velocity histograms for these galaxies in Fig- 
ure 5, and now focus on the new observations of And I, 
II, and III. 

In section 1 we noted that only one previous spec- 
troscopic study of stars in And I and And III exists 
(Guhathakurta, Reitzel, & Grebel 2000), and two such 
studies of And II exist, both being based on the same 
data set (Cote, Oke, & Cohen 1999; Cote et al. 1999). For 
And I and III, Guhathakurta, Reitzel, & Grebel (2000) 
were able to measure velocities for 29 and 7 dSph member 
stars, respectively, but with large velocity errors. Cote, 
Oke, & Cohen (1999) established accurate radial veloci- 
ties for seven member stars in And II and confirmed 35 
other members with very poor radial velocity accuracy, 
a v > 40 km s" 1 (Cote et al. 1999b; P. Cote 2006, pri- 
vate communication). The first aim of our survey is to 
significantly increase the numbers of stars with accurate 
radial velocities (a v < 10 km s _1 ) in these galaxies. In 
Figure 6 (bottom), we present radial velocity histograms 
for all objects for which we could measure a velocity as 
discussed earlier (grey histogram). For each of our And 

I, II, and III fields, there is a very clean detection of a 
kinematically cold spike representing stars that belong 
to the dSph. We also see a small population of field M31 
stars with a broad range of velocities. Considering that 
each of And I, II, and III are located well beyond the 
radius at which M31's halo dominates its inner spheroid 
(25 - 30 kpc, Guhathakurta et al. 2005; Kalirai et al. 
2006a), these field stars likely represent M31 halo mem- 
bers. Finally, there is a population of foreground Milky 
Way dwarf stars at small negative velocities. 

We establish membership of stars belonging to And I, 

II, and III using a three step process. This involves first 
isolating the sample of stars that have characteristics of 
giants at the distance of M31, second, eliminating any 
M31 substructure, and, third, cleaning the dSph pop- 
ulations of any M31 field stars in the halo. For the 
first cut, we can use photometric and spectroscopic diag- 
nostics to remove any Milky Way dwarf stars along the 
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Fig. 6. — Top — The distribution of likelihoods of individual stars being giants at the distance of M31 {Li > 0) or foreground Milky 
Way dwarfs {Li < 0) for each of our And I, II, and III sightlines. The likelihoods are constructed using five different photometric and 
spectroscopic diagnostics as discussed in § 3 and Gilbert ct al. (2006). We consider the giants to be those stars that are three times more 
likely to be M31 giants as compared to Milky Way dwarfs {Li > 0.5). Middle - CMDs for And I, II, and III illustrating our giant (filled 
points) and dwarf (open points) stars. The giants clearly form a tight RGB that follows the most metal-poor isochrones shown (see text). 
The filled points with redder colors that do not follow the locus of most giants are field M31 halo stars, as discussed in § 3. Bottom - 
Radial velocity histograms for each of our And I, II, and III sightlines show a dominant cold population at large negative velocities (the 
dSphs) as well as a smaller underlying broad distribution (M31 field halo RGB stars). Foreground Milky Way dwarfs are also seen at 
small velocities and are well separated from the dSph stars. The grey histogram shows all stars for which we measured a velocity whereas 
the black histogram illustrates our giant sample (dSph members and M31 field halo stars). As discussed in §3, the And I sightline is 
contaminated by M31 field substructure with a similar velocity to the dSph. This substructure is present in this Figure, and is removed to 
create our final member dSph RGB sample as illustrated in Figure 7. 



line of sight. The diagnostics used for this separation 
include radial velocity measurements, DD051 photom- 
etry, strength of the Nai doublet, position of the star 
in the CMD, and a comparison of the photometric and 
spectroscopic metallicities of each star. Details on the 
resolving power of each of these diagnostics to weed out 
Milky Way dwarf star contamination, as well as a care- 
ful step-by-step cookbook on how they are applied to the 
data set, are presented in Gilbert et al. (2006). As those 
methods were developed to separate out stars in M31's 
(kincmatically hot) field halo from Milky Way contami- 
nants, we note that we made two small changes in this 
application. First, the velocity and dispersion of the ex- 
pected dSph stars were set to the rough values of the 
peaks in Figure 6, and second the distance of each satel- 
lite (e.g., in calculating a photometric metallicity) was set 



to the known distances of And I (745 ± 24 kpc), And II 
(652 ± 18 kpc), and And III (749 ± 24 kpc) as reported 
in McConnachie et al. (2005). The uncertainty in the 
adopted distance of each satellite produces a very small 
effect on the overall classifications, given both the nature 
of the measurement and the relative insensitivity of the 
two affected diagnostics (i.e., metallicity comparison and 
position in the CMD) compared to others (e.g., radial 
velocity and DD051 photometry). The final likelihoods 
of a star being a giant or dwarf are illustrated in the top 
panel of Figure 6; values of Li that are positive are more 
likely to be giants and values with Li < are likely Milky 
Way dwarfs. Our final selection of giants are those stars 
that are three times more likely to be giants than Milky 
Way dwarfs (corresponding to L, > 0.5), as discussed in 
Gilbert et al. (2006) and Kalirai et al. (2006a). 
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Fig. 7. — The And I dSph is spatially located within identified 
halo substructure from the Giant Southern Stream (Ibata et al. 
2007; Gilbert et al. 2009). Left - The substructure is isolated from 
the dSph stars with a metallicity cut, where the more metal-rich 
material is associated with the stream debris. Right - As expected, 
these stars, shown as red open circles, are also located in the out- 
skirts of our spectroscopic masks which intersect at the center of 
And I. 



The distribution in radial velocity of giant stars in 
fields And I, II, and III are shown as the solid histogram 
in Figure 6 (bottom). Clearly, this separation has re- 
moved the obvious Milky Way dwarf stars at small nega- 
tive velocities, however, we note that a few dwarf out- 
liers are also removed at (or close to) the velocity of 
the dSph satellites. The distribution of the M31 giants 
(filled points) and Milky Way dwarfs (open points) are 
also shown on the CMD in Figure 6 (middle panels). 
The giants form a tight sequence in the CMD, stretch- 
ing from the tip of the RGB to our photometric limit at 
I ~ 22, whereas the dwarf stars occupy a broad range 
in color and luminosity depending on their spectral type 
and distance. The three solid red curves on these dia- 
grams are theoretical isochrones for an age of 12 Gyr and 
[a/Fe] = 0.0, with a metallicity of [Fe/H] = -2.31 (left), 
— 1.31, and —0.83, and —0.40 (right) from the Vanden- 
berg, Bergbusch, & Dowler (2006) models. The dashed 
curves are the same isochrones assuming [a/Fe] = +0.3. 
As we will quantify in § 5.3, each of these dSphs appear to 
contain a metal-poor population of stars with a modest 
dispersion. 

The second cut in our membership criteria only ap- 
plies to And I, which is situated at a spatial location in 
M31 that appears to contain substructure related to the 
Giant Southern Stream (Ibata et al. 2007). As Gilbert 
et al. (2009) show, the substructure in this field is located 
at roughly the same velocity as And I, but is much more 
metal-rich than the dSph stars. We can see the substruc- 
ture in the CMD of this dSph as an extention of metal- 
rich stars. This is illustrated more clearly in Figure 7, 
where we show the M31 giants (circles) and Milky Way 
dwarfs (crosses) in the And I field on an [Fe/H] vs u rac j 
plane (left panel - see § 5.3 for the details on how [Fe/H] 
is calculated). The reported photometric metallicity of 
the Milky Way dwarfs is meaningless since their distance 
is unknown. However, the metallicity of And I and the 
stream contamination is correct since these populations 
are both located at M31's distance. The diagram shows 
a bi-modality with a split at [Fe/H] ~ —0.9, where the 
bulk of the And I stars (filled points - see below) com- 
prise the more metal-poor, tight grouping of stars. To 



further illustrate that the metal-rich stars in this diagram 
are likely not members of And I, we look at their posi- 
tions relative to our two spectroscopic masks in the right 
panel. As expected, the stars associated with the stream 
contamination are found in the outskirts of the masks, 
far from And I, where the densities of the dSph and the 
field become more equal. We note that the King limiting 
radius of the dSph is 10', and most of these stars are lo- 
cated beyond this radius. Formally, we have selected the 
exact cut that separates the filled circles from the open 
circles arbitrarily at the apparent break point in Figure 7 
(left). For our sample, this translates to a hard cut of 
[Fe/H] < —0.92 for the And I members, although we rec- 
ognize that the exact separation may be slightly shifted 
from this value and depends on the relative age and dis- 
tance difference between the two populations. The hard 
cut may also lead to the truncation of any metal-rich tail 
in the dSphs metallicity distribution function. 

The final cut in our membership criteria eliminates 
any likely M31 field halo giants not associated with each 
dSph. The limiting radii of And I is 10.4', of And II is 
22.0', and of And III is 7.2' (McConnachie & Irwin 2006), 
and so we eliminate three stars in And I and eight stars 
in And III that are located beyond these limits. These 
eliminations work to actually remove stars that can be 
easily identified as outliers on the velocity histograms. 
For example, in And I one of the stars eliminated has 
a radial velocity of —440.0 km s _1 , and in And III one 
of the stars is moving at u ra d = —431.6 km s _1 . Both 
of these stars have velocities that are 5 - 10a deviant 
from the final observed dispersion. Additionally, the cut 
at the limiting radius in And III eliminates the two gi- 
ants with V-I > 2 in the CMD in Figure 6. The pho- 
tometric metallicity of these two stars is calculated to 
be [Fe/H] > —0.25, well over 10 times more metal-rich 
than the remaining dSph stars (see § 5.1). The resulting 
velocity distribution still includes M31 field halo giants 
that are spatially superimposed on the dSph population. 
To eliminate these last few stars, we measure the mean 
and sigma of the resulting velocity distribution and iter- 
atively eliminate stars that are more than 3-c from the 
mean. This iteration removes just two stars from And I, 
zero stars from And II, and three stars from And Ill's 
sightline. The eliminated stars were again not centrally 
located in the dSph suggesting they are likely field halo 
giants. The final selected members of these three dSphs 
include 80 stars in And I, 95 stars in And II, and 43 stars 
in And III. 

We stress that our primary goal from the selection pro- 
cess described above is to construct as secure a sample 
as possible of confirmed RGB members of these three 
satellites. Given the cuts used, it is possible (although 
unlikely) that we have in fact thrown a member or two 
out of the sample, however our sample sizes are large 
enough that the exclusion of a few stars will not impact 
the analysis that follows. The filled histograms in Fig- 
ure 8 (bottom) illustrate the velocities for the stars in 
the final sample, and the filled points in the top panel 
of this Figure show their distribution in the CMD. The 
open histograms (bottom) and open points (top) in Fig- 
ure 8 represent the eliminated giants from the limiting 
radius and c-clipping cuts discussed above. 
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Fig. 8. — A closer view of the CMD (top) and radial velocity histograms (bottom) for our final And I (80 stars), And II (95 stars), and 
And III (43 stars) samples of RGB members (filled points and histograms). Open points and histograms show the M31 field halo giant 
stars that are eliminated from the final sample as discussed in § 3. The three sets of isochrones shown in each panel of the CMDs are taken 
from Vandenberg, Bergbusch, & Dowler (2006), for an age of 12 Gyr and metallicities of [Fe/H] = —2.31 (left), —1.31 (middle), and —0.83 
(right), and have been shifted to each satellites' distance modulus. The models for [a/Fe] = 0.0 are shown as solid curves and those for 
[a/Fe] = +0.3 are shown as dashed curves. Clearly, all three of these satellites are metal-poor systems. 



4. MEAN VELOCITIES, INTRINSIC VELOCITY 
DISPERSIONS, AND M/L RATIOS 

We calculate the mean radial velocity and intrinsic dis- 
persion for each of the six dSphs using the maximum- 
likelihood method described by Walker et al. (2006). 
This method assumes that the observed dispersion is the 
sum of the true intrinsic dispersion and the velocity er- 
rors discussed at length in §2.4. For And I, we find w ra( j 
= -375.8 ± 1.4 km s" 1 and a v = 10.6 ± 1.1 km s -1 , 
for And II we find D ra( j = —193.6 ± 1.0 km s _1 and 
a v = 7.3 ± 0.8 km s -1 , for And III we find u ra( j = 
-345.6 ± 1.8 km s" 1 and a v = 4.7 ± 1.8 km s _1 , for 
And VII we find v ra d = —309.4 ± 2.3 km s _1 and a v — 
9.7 ± 1.6 km s" 1 , for And X we find w rad = -168.3 ± 

1.2 km s^ 1 and a v = 3.9 ± 1.2 km s _1 , and for And XIV 
we find t> ra d = —481.0 ± 1.2 km s _1 and a v = 5.4 ± 

1.3 km s . The mean radial velocities of these six galax- 
ies are consistent with M31 membership, since M31 has 
a systemic velocity of —300 km s _1 . The one exception 
to this may be And XIV, which has both a large negative 
radial velocity and a large distance from M31 (both pro- 
jected and along the line of sight). As discussed in Ma- 
jewski et al. (2007), the galaxy may therefore be falling 
into the Local Group for the first time. The velocities for 
And I and II are inconsistent with the measurements by 
previous studies within the mutual 1-a errorbars estab- 
lished by this and the previous studies (Guhathakurta, 
Reitzel, & Grebel 2000; Cote et al. 1999). As we noted 
earlier, these previous results are only based on a handful 



of stars and therefore our values are much more secure. 
For And III, Guhathakurta, Reitzel, & Grebel (2000) 
found w ra d = —352.3 ± 13.6 km s -1 based on 9 RGB 
members, which is consistent with our much more pre- 
cise result. 

The central velocity dispersion of the six galaxies in 
this study varies by more than a factor of 2.5. The only 
previous measurement to compare our results with is the 
high resolution study of Cote et al. (1999), which estab- 
lished a v = 9.3 ± 2.7 km s _1 for And II based on seven 
RGB stars. Our result is consistent within the larger 
error bar determined by Cote et al. (1999) for their mea- 
surement, however we find that the velocity dispersion 
is lower. We can use our new velocity dispersions to 
put first order constraints on the total mass (stellar and 
dark matter) of each satellite using the method described 
by Illingworth (1976), M = l&7*Pr c a 2 v where j3 = 8 for 
dSphs (Mateo 1998). This formalism was initially con- 
structed for the analysis of line of sight velocities in glob- 
ular clusters and therefore has several assumptions built 
into it, all of which may not be true. For example, we 
are assuming the galaxies are spherical, are in dynamical 
equilibrium, and have an isotropic velocity dispersion. 
We are also assuming that the stellar distribution fol- 
lows a King profile (which is observed) and traces the 
dark matter. We return to a discussion on how more ac- 
curate dynamical mass measurements can be established 
for M31 dSphs below. 

The core radius of four of the six dSphs is known from 
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the structural study of McConnachie & Irwin (2006), who 
find r c = 580 ± 60 pc for And I, r c = 990 ± 40 pc for 
And II, r c = 290 ± 40 pc for And III, and r c = 450 ± 
20 pc for And VII (geometric means). The measured core 
radius, and an assumed 10% uncertainty, of And X is r c 
= 270 ± 30 pc (Zucker et al. 2007) and that of And XIV is 
much larger at r c — 730 ± 70 pc (Majewski et al. 2007). 
The total masses of these satellites are therefore M = 
(8.7 ± 1.6) x 10 7 M Q for And I, M = (7.0 ± 1.1) x 10 7 
M for And II, M = (8.6 ± 4.8) x 10 6 M for And III, 
M = (5.7 ± 1.3) x 10 7 M Q for And VII, M = (5.5 ± 2.5) 
x 10 6 M Q for And X, and M = (2.9 ± 1.0) x 10 7 M G for 
And XIV. At the luminous end, for M v = -11 to -13, 
the masses of And I, II, and VII are similar to Milky Way 
dSphs (e.g., Fornax and Leo I), calculated in the same 
way. At lower luminosities, the masses of And III and X 
appear to be lower than systems in the Milky Way with 
similar brightness, such as Carina, Draco, and Sextans. 

We can calculate the mass-to-light ratios in Solar units 
of And I, II, III, VII, X, and XIV by combining the 
masses determined above with the measured luminosities 
of each satellite (see McConnachie & Irwin 2006; Zucker 
et al. 2007; Majewski et al. 2007). This gives M/L v = 
19 ± 4 for And I, M/L v = 7.5 ± 1.8 for And II, M/L v 
= 8.3 ± 5.2 for And III, M/L v = 3.2 ± 1.2 And VII, 
M/Ly = 37 ± 24 for And X, and M/L v = 160 ± 95 for 
And XIV. The mass-todight ratios indicate that these 
satellites are dark matter dominated, as expected, how- 
ever both And II and III are at the low end of the range 
of M/L ratios for dSphs. A more in depth comparison of 
the results in this section, as well as their implications, 
will be presented in § 6.4. 

In Table 2 we summarize several of the fundamental 
properties of each of the six satellites discussed above. 
This includes their luminosities, projected distance from 
M31, mean radial velocities, intrinsic velocity disper- 
sions, metallicities, and metallicity dispersions (see § 5 
for half-light radii, total masses, and M/L ratios). Fur- 
ther details on the individual galaxies are also available 
in the Guhathakurta, Reitzel, & Grebel (2000), Kalirai 
et al. (2009), and Majewski et al. (2007) studies. 

4.1. More Accurate Masses and Mass 
Profiles of M31 's dSphs 

Both Wolf et al. (2009) and Walker et al. (2009a) have 
recently compared the masses of Milky Way dSphs cal- 
culated using the method above to those derived using 
new mass estimators for dispersion-supported systems 
that result from a manipulation of the Jeans equation. 
The new mass estimators are shown to place tight con- 
straints on the system mass within the half light radius, 
and weaker constraints at larger radii. In the Walker 
et al. (2009a) analysis, the mass at the half light radius 
using the new formalism is found to be less than the 
Illingworth (1976) (total mass) approximation by ~50% 
(systematically). Wolf et al. (2009) go one step further 
and compare the estimated total masses of dSphs (see 
below) to the Illingworth (1976) relation, and find that 
the latter underpredicts both the total mass of dSphs 
and the associated uncertainties in the mass. The reason 
for this is related to the assumed distribution of the total 
mass, where the Illingworth (1976) approximation forces 
the mass distribution to truncate near the stellar extent 
of the galaxy. However, over a range of almost two or- 



ders of magnitude in mass, the offset in mass is nearly 
constant (see their Figure CI in the Appendix). 

The baseline for the comparisons above that demon- 
strates the mass at the half light radius is well con- 
strained in these new derivations comes from the full 
modeling of individual radial velocities to yield mass pro- 
files for Milky Way dSphs. For example, Strigari et al. 
(2008b) assume the radial velocities are related to the 
overall mass distribution through the Jeans equation, 
and that the dark matter follows a five parameter density 
profile. They also allow the velocity anisotropy to vary in 
the modeling, however, the systems are still assumed to 
be spherically symmetric and in dynamical equilibrium 
(see also Walker et al. 2009a). Given this recent work, 
we calculate the mass at the half light radius for each of 
the M31's dSphs using the new formalism presented in 
Wolf et al. (2009), M 1/2 = 3G" 1 (<^>r 1/2 , where r 1/2 is 
the three dimensional deprojected half light radius (see 
Table 2). This calculation yields M 1/2 = (7.0 ± 1.2) x 
10 7 M Q for And I, M 1/2 = (6.1 ± 1.0) x 10 7 M for 
And II, M 1/2 = (9.6 ± 5.4) x 10 6 M & for And III, M 1/2 
= (6.9 ± 1.6) x 10 7 M & for And VII, M 1/2 = (4.7 ± 2.0) 
x 10 6 M Q for And X, and M 1/2 = (9.2 ± 4.4) x 10 6 M 
for And XIV. These results are summarized in Table 2 
and discussed further in § 6.4. Additional analysis of the 
mass profiles of these satellites in comparison to Milky 
Way dSphs will also be presented in Wolf et al. (2010, in 
prep.), 

5. CHEMICAL ABUNDANCES 

Da Costa et al. (1996; 2000; 2002) measured the abun- 
dances of And I, II, and III from #ST/WFPC2 photo- 
metric data and found that all three dSphs are metal- 
poor, ([Fe/H] < —1.45), but with large variations in the 
internal abundance spread. For And I, they found a very 
high dispersion of o [fc/h] — 0.60, whereas for And III, 
they find that cpe/H] = 0.12. For the more luminous 
satellite And II, Da Costa et al. find (T[f c /h] = 0.36. 
These large variations suggest that the M31 dSphs may 
have experienced quite different evolutionary histories. 
For example, a galaxy such as And III would have had a 
tough time retaining its enrichment products whereas a 
system such as And I may have experienced cither multi- 
ple epochs (or an extended epoch) of star formation. Da 
Costa et al. stress the need for an independent spectro- 
scopic study of RGB stars in these satellites. As we've 
already demonstrated in § 3, the And I sightline contains 
bright red giants that belong to the debris of the Giant 
Southern Stream. These stars are much more metal-rich 
than the dSph population, and if not removed, will arti- 
ficially inflate the measured abundance spread. 

As in our previous studies of M31's stellar halo (e.g., 
Kalirai et al. 2006a), we determine the chemical abun- 
dance and abundance spread of stars in And I, II, and 
III using two independent methods. Both methods are 
used on the restricted sample of confirmed RGB mem- 
ber stars in each satellite, as determined in §3. The 
first method is based on a comparison of the location of 
the dSph stars on the (I, V—I) CMD to a new set of 
theoretical isochrones. The second method is based on 
a spectroscopic measurement of the equivalent width of 
the Can triplet (A <~ 8500 A). The advantages of this 
two-pronged photometric and spectroscopic approach is 
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TABLE 2 



Property And I And II And III And VII And X And XIV 

M Y -11.8 ± 0.1 -12.6 ± 0.2 -10.2 ± 0.3 -13.3 ± 0.3 -8.1 ± 0.5 -8.3 ± 0.5 

■Rproj (kpc) 44.9 145.6 68.2 214.5 75.5 162.5 

*) rad (km s- 1 ) -375.8 ± 1.4 -193.6 ± 1.0 -345.6 ± 1.8 -309.4 ± 2.3 -163.8 ± 1.2 -481.0 ± 2.0 

Uv (km s- 1 ) 10.6 ± 1.1 7.3 ± 0.8 4.7 ± 1.8 9.7 ± 1.6 3.9 ± 1.2 5.4 ± 1.3 

[Fe/H] -1.45 ± 0.04 -1.64 ± 0.04 -1.78 ± 0.04 -1.40 ± 0.30 -1.93 ± 0.11 -2.26 ± 0.05 

a [Fe/H] 0.37 0.34 0.27 •■■ 0.48 

r-c 1 580 ± 60 990 ± 40 290 ± 40 450 ± 20 271 ± 27 734 ± 73 

2D R e (pc) 2 682 ± 57 1248 ± 40 482 ± 58 791 ± 45 339 ±6 413 ± 41 

3D n/ 2 (pc) 3 900 ± 75 1659 ± 53 638 ± 77 1050 ± 60 448 ± 8 461 ± 155 

M (M ) 4 (8.7 ± 1.6) x 10 7 (7.0 ± 1.1) x 10 7 (8.6 ± 4.8) x 10 6 (5.7 ± 1.3) x 10 7 (5.5 ± 2.5) x 10 6 (2.9 ± 1.0) x 10 7 

M/L(M Q /L Q ) 4 19 ± 4 7.5 ± 1.8 8.3 ± 5.2 3.2 ± 1.2 37 ± 24 160 ± 95 

M 1/2 (M ) 5 (7.0 ± 1.2) x 10 7 (6.1 ± 1.0) x 10 7 (9.6 ± 5.4) x 10 6 (6.9 ± 1.6) x 10 7 (4.7 ± 2.0) x 10 6 (9.2 ± 4.4) x 10 6 

M 1/2 /L 1/2 (Mq/Lq) 5 31 ± 6 13 ± 3 19 ± 12 7.7 ± 2.8 63 ± 40 102 ± 71 

No. Stars 80 95 43 18 22 38 



1 r c is the core radius, the uncertainty is taken to be 10% for 
And X and And XIV. 

2 R e is the two dimensional elliptical half light radius. 

3r l/2 is tne three dimensional deprojectcd half light radius. 

4 Thcse masses represent the total system mass, as measured using 
the Illingworth (1976) formalism. 

5 The masses at the half light radius are calculated using M 1 / 2 = 
3G- 1 (al)r 1/2 (Wolf ct al. 2009). 



that (1) a comparison of two independent measures can 
shed light on the existence of a systematic bias in either 
method, (2) the radial velocity measurements from the 
spectroscopic data can be used to ensure that the photo- 
metric and spectroscopic abundance distribution consists 
only of true dSph RGB stars (i.e., they are uncontami- 
nated by Milky Way dwarfs) , and (3) the spectra for like 
stars can be coadded to yield information on detailed 
chemical abundances (e.g., [a/Fe]). Below, we discuss 
each of these two methods in turn to derive the abun- 
dances of stars in And I, II, and III. 

5.1. Photometric Metallicity Determination 

The positions of stars along the RGB in a stellar pop- 
ulation can be used to determine their chemical abun- 
dances, provided their age is known. In an (/, V— I) 
CMD, the shape of the RGB is such that metal-rich stars 
become increasingly redder relative to metal-poor stars 
as their luminosity increases. The stars that we targeted 
in And I, II, and III are among the brightest RGB stars 
in these dSphs, and therefore we take full advantage of 
this color sensitivity. 

The spectroscopically selected final sample of member 
stars in And I, II, and III are shown as filled circles in 
the top panel of Figure 8. As discussed above, the red 
curves represent theoretical stellar isochrones from the 
new models of Vandcnbcrg, Bergbusch, & Dowler (2006), 
for three different metallicities, [Fe/H] = —2.31 (left), 
— 1.31 (middle), and —0.83 (right). The solid curves 
are models with no alpha-enhancement and the dashed 
curves assume [a/Fe] = +0.3. The models have been 
shifted to the known distance to each satellite, as mea- 
sured from the tip of the red giant branch. We measure 
the metallicity distribution function for each satellite by 
interpolating the magnitude and color of each member 
star within a larger grid of similar isochrones, includ- 
ing a dozen models within the range [Fe/H] = —2.31 to 



—0.83. For the few stars in each satellite that are located 
outside the bounds of these models, we extrapolate the 
metallicity measurements. The errors on these measure- 
ments therefore depend both on the distance uncertainty 
to the dSph and the error in the photometry. 

The photometric metallicity measurements of And I, 
II, and III assume a particular set of models, in this 
case the new Vandenberg, Bergbusch, & Dowler (2006) 
isochrones. We can also recalculate the metallicities as- 
suming a different isochrone set. For example, our re- 
sults agree with the metallicities that would be derived 
from either the Padova models (Girardi et al. 2002) or 
the Yale-Yonsei models (Y 2 ; Demarque et al. 2004) at 
the 0.15 dex level. The metallicities also depend on the 
assumed age of the stellar population. For intermediate 
to old ages, there is little color sensitivity of an RGB 
star with age. For example, a shift in the entire age 
of the dSph from 12 Gyr (our adopted age) to 6 Gyr 
would translate to a +0.2 dex offset in [Fe/H] p hot- Al- 
though there currently exist no direct age measurements 
for And I, II, or III from resolved main-sequence turnoff 
fitting, several hints suggest that the bulk of the stellar 
populations in each of these three dSphs are old. For ex- 
ample, the CMDs in Da Costa et al. are similar to those 
of Galactic globular clusters, which are all old. Specifi- 
cally, the galaxies show an absence of upper asymptotic 
giant branch stars, contain blue horizontal branch stars, 
and contain RR Lyrae stars. 

5.2. Spectroscopic Metallicity Determination 

The spectroscopic metallicity of a star, which we des- 
ignate [Fe/H] spoc , can be derived from the equivalent 
widths of the three Can absorption lines. The inde- 
pendent measurements of the three lines are combined 
to produce a reduced equivalent width as described in 
Rutledge, Hesser, & Stetson (1997). This reduced equiv- 
alent width is then empirically calibrated based on ob- 
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servations of Galactic globular cluster RGB stars to yield 
[Fe/H] spcc (Rutledge et al. 1997). The correction for sur- 
face gravity is derived from the luminosity of the HB 
in the HST imaging study of each dSph, which we take 
from the Da Costa et al. work. We measure this from 
their CMDs to be V HB = 25.20 for And I, V KB = 24.85 
for And II, and Vhb = 24.95 for And III. 

The uncertainty in [Fe/H] spoc from our low S/N spec- 
tra results in a larger scatter than the [Fe/H] phot mea- 
surements. Still, we can use this measurement to provide 
an independent check on our overall metallicity measure- 
ments. In Figure 9, we compare the two metallicity mea- 
surements to one another. We have shown all individ- 
ual data points, measured using both techniques, as well 
as binned averages of the spectroscopic measurements 
(larger filled circles with error bars). These were calcu- 
lated using a minimum bin size of 0.2 dex in [Fe/H] p hot 
while ensuring >15 stars in each bin, and requiring the 
spectra of the individual stars to have at least S/N > 
3 (see below). Out of our initial sample of 218 member 
stars in And I, II, and III, this S/N cut leaves 158 objects. 
We reduce this sample by three more stars to eliminate 
objects for which the [Fe/H] spoc values were unrealisti- 
cally high (+0.71, +2.05, and +2.61). There is a nice 
overall agreement between the resulting mean values of 
the two measurements over a range that includes most of 
the data points, [Fe/H] p hot = —2.2 to —1.1. The dashed 
line illustrates the 1:1 relation. For our most metal-poor 
bin, [Fe/H] p hot < —2.2, the spectroscopic metallicities 
are systematically more metal-rich than the photometric 
values. For stars with these low metallicities, the Call 
lines are very shallow and the spectroscopic metallicities 
using this method are likely in error. In fact, Kirby et al. 
(2008) show that the spectroscopic metallicities of sim- 
ilar stars in Galactic dSphs have been overestimated in 
past studies. 

We also illustrate how the independent photometric 
and spectroscopic metallicity measurements compare for 
various S/N cuts. The crosses in Figure 9 illustrate those 
objects with S/N < 3, the open circles represent stars 
with 3 < S/N < 5, and the filled points are those with 
S/N > 5. The higher S/N points tend to be clustered 
closer to the 1:1 line as compared to the measurements 
from poorer quality spectra, indicating that much of the 
vertical scatter in the diagram results from poor charac- 
terization of the Can triplet lines in these spectra. 

To summarize, we find a nice overall agreement be- 
tween our photometric and spectroscopic metallicity 
measurements for these satellites, indicating that any 
systematic biasses are likely small. For example, if we 
assume the dSphs are significantly younger than 12 Gyr, 
the resulting comparison of metallicities leads to a re- 
lation offset than the 1:1 line. We illustrate this, for 
an age of 6 Gyr, as the green curve in Figure 9. We 
also investigate the effects of different a-enhancements. 
The photometric metallicities shown in Figure 9 assume 
these three dSphs are not enhanced in a-elements rela- 
tive to the Sun (e.g., we used the isochrones with [a/Fe] 
= 0.0 for these calculations). Stars in the Milky Way's 
field halo, and globular clusters, are known to be a- 
enhanced ([a/Fe] = +0.3) and it is generally believed 
that this is a result of the stars in the halos of galax- 
ies forming in early "bursts". Milky Way dSphs also 
show some a-enhancement, although not as high as that 
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Fig. 9. — The [Fe/H] p hot measurements are compared to the 
independently measured [Fe/H] sp cc values, over a 1.5 dex spread 
in [Fe/HJphot- The dashed line illustrates equality between the two 
measurements, and, in general, our data points follow this trend 
nicely for [Fe/H] p hot > —2.2. The scatter among the spectroscopi- 
cally measured abundances is large, as expected given the low S/N 
of our spectra. This is illustrated by plotting the symbols as a func- 
tion of the spectral S/N, where crosses represent all stars with a 
spectral S/N < 3, open circles represent stars with 3 < S/N < 
5, and filled circles represent stars with S/N > 5. The scatter in 
[Fe/H] spcc is clearly reduced among the better quality data. The 
larger data points with error bars are binned averages of the S/N 
> 3 data, where the photometric metallicities have been calculated 
using the Vandenberg, Bergbusch, & Dowler (2006) isochrones with 
an age of 12 Gyr and [a/Fe] = 0.0. The two solid blue and green 
lines show the same comparison for different assumptions on age 
and [a/Fe], as discussed in §5.2. 



measured in globular clusters. For example, Shetrone, 
Cote, & Sargent (2001) find +0.02 < [a/Fe] < +0.13 
for the three dSphs Draco, Sextans, and Ursa Minor. A 
mild a-enhancement, such as that seen in these Milky 
Way galaxies, produces a very small affect on our metal- 
licity measurements. However, we note that if And I, 
II, and III are enhanced in a-elements at a level similar 
to Galactic globular clusters (e.g., [a/Fe] = +0.3), then 
our overall photometric metallicities would become more 
metal-poor by ~0.1 - 0.2 dex. This comparison is shown 
as a blue curve in Figure 9. 

Our two independent metallicity measurements gener- 
ally agree with one another, and we have shown that the 
scatter in the spectroscopic metallicities results largely 
from poor S/N. Therefore, we will adopt the photomet- 
ric metallicities for the subsequent analysis. We stress 
that the reported abundances are derived from a secure 
sample of member stars in each satellite. 

5.3. Abundance Distributions 

The photometric metallicity distribution functions 
(MDFs) for the confirmed RGB stars in each of And I, 
II, and III in displayed in Figure 10. The top panels 
show discrete histograms and the bottom panels show 
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Fig. 10. — The discrete photometric metallicity distribution functions (MDF) for And I, II, and III (top) and the corresponding cumulative 
distributions (bottom), based on our confirmed sample of member stars in each dSph. The dashed line is held fixed at [Fc/HJphot = —1.75 
in each panel as a guide. We find that And III is the most metal-poor dSph of the three, and that And I is the most metal-rich. The 
metallicity dispersion in all three dSphs is similar, <T[p e /H] = 0.3 - 0.4. 



the cumulative distributions. The dashed line marks a 
fixed guide at [Fe/H] = —1.75 in both panels. As was 
apparent from the CMDs, we confirm that all three of 
these dSphs are in fact metal-poor galaxies, with And I 
being slightly more metal-rich relative to And II and III. 
We find that the internal abundance spread of all three 
galaxies is quite similar, differing by less than 0.1 dex. 
Formally, we measure the mean metallicity, error in the 
mean, and dispersion to be [Fe/H] p hot = —1-45 ± 0.04 
(qFe/H] = 0.37) for And I, [Fe/H] phot = -1.64 ± 0.04 
(ff[ Fc /H] = 0.34) for And II, and [Fe/H] phot = -1.78 ± 
0.04 [Fc/H ] = 0.27) for And III. 

Abundance analysis of And VII, X, and XIV were 
presented in Grebel & Guhathakurta (1999), Kalirai 
et al. (2009), and Majewski et al. (2007). Summariz- 
ing, the high resolution spectroscopic data for And VII 
in Guhathakurta, Reitzel, & Grebel (2000) did not yield 
an abundance measurement due to the low S/N of the 
18 confirmed giants at the Can triplet. Therefore, we 
adopt And VIFs metallicity from the Keck/LRIS pho- 
tometric study of Grebel & Guhathakurta (1999), who 
find [Fe/H]phot = —1-4 ± 0.3 (the dispersion in metallic- 
ity was not measured). For And X, Kalirai et al. (2009) 
measured both a spectroscopic and photometric metal- 
licity and found these to be in nice agreement. They 
find that And X is a metal-poor galaxy with [Fe/H] p hot 
= —1.93 ± 0.11, and exhibits a large metallicity disper- 
sion of ct[ Fo / H ] = 0.48. For a distance modulus of 24.7, 
the the photometric metallicity of And XIV is very low, 



[Fe/H] = —2.26 ± 0.05 (there is no metallicity disper- 
sion measurement). These results are also summarized 
in Table 2. 

For And I, II, and III, we can compare our metallic- 
ity results to the Da Costa et al. (1996; 2000; 2002) 
#ST/WFPC2 study. Da Costa et al. transformed their 
photometry from the native HST filters to Johnson- 
Cousins, and then compared the color of the RGB stars 
to a set of giant branches for Galactic globular clusters 
(i.e., alpha-enhanced populations). As we noted earlier, 
their study of And I is likely affected by the metal-rich 
stream debris that we eliminated from this sightline in 
§ 3. Interestingly, our mean metallicity is still identical to 
what they found, [Fe/H] p hot = —1-45 ± 0.20, suggesting 
that their value was in fact underestimated. As expected, 
our dispersion is much smaller than their study, CT[Fe/H] 
= 0.60. Our sample only includes stars with [Fe/H] < 
—0.92 given the detected substructure in this field, and, 
therefore, we would have eliminated the presence of any 
metal-rich tail in this dSph. For And II, we find a very 
similar abundance spread in the galaxy compared to the 
results of Da Costa et al. (ct[f c /h] = 0.36), however our 
mean metallicity is more metal-poor than their study by 
^0.15 dex. For And III, our mean metallicity is consis- 
tent with the Da Costa et al. study within their larger 
error bar ([Fe/H] p hot = —1.88 ±0.11), however our in- 
ternal metallicity dispersion is more than 2 x larger than 
their measurement (a = 0.12). 

Recently, McConnachic, Arimoto, & Irwin (2007) im- 



16 



Kalirai et al. 



aged And II using the wide-field Subaru Suprime-Cam 
instrument to a photometric depth below the horizontal 
branch. Their analysis suggests that the dSph contains 
two distinct components. For the dominant extended 
component, they find an old population with [Fe/H] = 
— 1.5 and cr[Fc/H] = 0.28. For the more concentrated in- 
ner region of the dSph, they find an intermediate aged 
(7 - 10 Gyr) and metal-rich ([Fe/H] = —1.2 and cr[ Fe /H] 
= 0.40) population. Confirming this picture with our 
present spectroscopic data is difficult given the limited 
spatial coverage (e.g., lack of many stars near the cen- 
ter of the galaxy), however, we have now obtained five 
times as many spectroscopic measurements over the face 
of this galaxy. These data, consisting of ^500 individual 
radial velocity measurements of And II members, will 
present a much cleaner view of the nature of And IPs 
stellar populations. 

Overall, our results suggest that the chemical abun- 
dances and abundance spread (where available) of the 
brighter four M31 satellites in our sample are quite sim- 
ilar, and that there is no evidence suggesting the evo- 
lutionary histories of these galaxies were significantly 
different. The two faintest satellites in this clean sam- 
ple, And X and XIV, are found to be more metal-poor 
than the brighter satellites. Of course, a full analysis 
will require deeper photometric observations of the main- 
sequence turnoff morphology in each system, from which 
a full star formation history can be derived in conjunc- 
tion with these metallicity measurements. We also note 
that the metallicities of the M31 halo dSphs is similar 
to the mean metallicity of M31's stellar halo, determined 
by Kalirai et al. (2006a) and Chapman et al. (2006) to 
be [Fe/H] ~ -1.5. 

6. GLOBAL PROPERTIES: 
COMPARING MILKY WAY AND M31 DSPHS 

6.1. An Inventory of the Milky Way and M31 dSphs 

As we introduced earlier, wide field imaging surveys 
have recently uncovered many new dSph galaxies orbit- 
ing the Milky Way. As rapidly as these systems are 
being discovered, different groups have targeted individ- 
ual giant and dwarf stars with multiobject spectrographs 
to characterize the radial velocities, velocity dispersions, 
and spectroscopic abundances of these systems. For com- 
parison to our M31 sample, we begin by considering only 
those Milky Way dSphs in which the intrinsic velocity 
dispersion has been resolved. 

Our baseline catalogue of properties for the classical 
Milky Way dSphs is drawn from Irwin & Hatzidimitriou 
(1995) and the review by Mateo (1998), updated to re- 
flect recent photometric and spectroscopic analysis on 
a galaxy by galaxy basis where available. Our sample 
includes the Milky Way dSphs Carina, Draco, Fornax, 
Leo I, Leo II, Sculptor, Sextans, and Ursa Minor, as well 
as the more distant Local Group dSphs Tucana and Ce- 
tus (McConnachie & Irwin 2006; Lewis et al. 2007). For 
the core radii and central velocity dispersions, we adopt 
values from Wolf et al. (2009) who present a summary 
of the most updated results in their Table 1 (see full 
references in their study, including Walker, Mateo, & 
Olszewski 2009). The [Fc/H] of most of these galaxies 
is updated from the Mateo (1998) results based on new 
multiobject studies of larger numbers of stars. For exam- 
ple, for Carina, Fornax, Sculptor, and Sextans we use the 



latest DART team results (Helmi et al. 2006), for Leo I 
and Leo II we use the Koch et al. (2007a) and Koch et al. 
(2007b) studies, and for Draco and Ursa Minor we adopt 
the values in Winnick (2003). If it is not directly mea- 
sured in these studies, we have assumed an uncertainty 
of 10% in the half light radius of each galaxy. For the 
newly discovered lower luminosity Milky Way dSphs, our 
sample includes Bootes I, Bootes II, Canes Venatici I, 
Canes Venatici II, Coma Berenices, Hercules, Leo IV, 
Leo T, Segue I, Ursa Major I, and Ursa Major II. The 
baseline properties, such as luminosity and half-light ra- 
dius, come from the summary table in Wolf et al. (2009) 
(see references therein and also the structural analysis in 
Martin, de Jong, & Rix 2008), velocity dispersions come 
from the Keck/DEIMOS spectroscopic study by Simon 
& Geha (2007), and abundances come from Kirby et al. 

(2008) . The core radii of the dSphs is again taken from 
the King model fits in Strigari et al. (2008b), with a 10% 
assumed error. We note that the core radius of Bootes I 
is not known and so we have made an approximation 
such that r c ~ R e . This can lead to a ^25% uncertainty 
(Martin et al. 2007), which has been factored into the 
calculation of dynamical masses (see below). 

For the 19 M31 dSphs, plus And XVIII which appears 
to be a background object, we adopt parameters primar- 
ily from the discovery papers summarized earlier as well 
as the analysis of McConnachie & Irwin (2006) for the 
brightest satellites. The half light radii in McConnachie 
& Irwin (2006) are computed as geometric means, and so 
we have adjusted these to an elliptical half-light radius 
(R e ) by factoring in the measured cllipticity. The metal- 
licity of And V is taken from Armandroff, Davies, & Ja- 
coby (1998) and for And VI and VII is taken from Grebel 
& Guhathakurta (1999). For the six satellites we've dis- 
cussed in this paper, we adopt our new results based 
on the spectroscopically clean samples of stars. Again, 
a 10% error in the core radius is assumed (i.e., for the 
calculation of dynamical masses) if it has not been mea- 
sured in the earlier papers. As discussed in § 1, the lim- 
ited kinematic data presented for And IX and And XII 
(Chapman et al. 2005, 2007; Collins et al. 2009), And XI 
and XIII (Collins et al. 2009), and And XV and XVI 
(Letarte et al. 2009) have large uncertainties and are in- 
conclusive in yielding accurate dynamical masses, and so 
we ignore them in the present analysis. 

We recalculate the dynamical mass of each of the Milky 
Way dSphs summarized above using the same (simpli- 
fied) technique described in § 4, with updated veloc- 
ity dispersion and core radius measurements from these 
studies. More accurate masses of these Milky Way satel- 
lites at the half light radii are calculated in Wolf et al. 

(2009) . A summary of all of these current properties, 
for both the Milky Way dSphs with internal velocity dis- 
persion measurements, and for all of the M31 dSphs, is 
given in Table 3. 

6.2. Chemical Abundance Trends 

Table 3 illustrates that the Milky Way and M31 each 
host ~20 dSphs, and the bulk of these satellites have 
measured abundances and half light radii. The intrinsic 
luminosity range of the M31 dSphs is of course smaller 
than the Milky Way satellites, however the two distribu- 
tions overlap over a factor of 1000 in luminosity. There 
exists a known luminosity-metallicity relation for Local 
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Fig. 11. — The properties of M31 dSphs are compared to Milky Way dSphs of the same luminosity (see §6 for references). The Milky 
Way dSphs are illustrated with as grey points whereas the M31 points are darker. The six M31 galaxies for which we have resolved intrinsic 
velocity dispersion measurements from the present study are dark encircled points. Finally, the three distant Local Group dSphs, Cetus, 
Tucana, and And XVIII are shown with dotted points and errorbars. All results are summarized in Table 3. The left panel illustrates 
the luminosity-metallicity relation, where we find excellent overlap between the M31 and Milky Way satellites, suggesting these systems 
have had similar chemical evolution histories. Four of the faintest M31 dSphs are found to be signicantly more metal-rich than expected 
from the trend, and therefore may be systems that have experienced tidal stripping (i.e., they are underluminous for their metallicity) . 
The error bars in this plot represent the dispersion in metallicity in each galaxy. For several dSphs lacking such measurements (mostly in 
M31), the dispersion has been arbitrarily set to 0.30 dcx (sec Table 3). The middle-left panel illustrates the luminosity-size relation, where 
the size is represented by the 2D elliptical half-light radius. Although the brightest M31 satellites are larger than Milky Way dSphs of the 
same luminosity (see also McConnachie & Irwin 2006), we find that lower luminosity systems have similar sizes and therefore dynamical 
effects are unlikely to have truncated the sizes of these systems in a systematic way. The masses and mass-to-light ratios of the M31 dSphs, 
as calculated using the simplified technique described in § 4, are illustrated in the two right panels. The overall trend of decreasing total 
mass-to-light ratio with increasing luminosity is similar for the dSphs of each galaxy, although two of the faintest M31 dSphs in our sample 
appear to be less massive than their Milky Way counterparts. Combined with their similar (or larger) sizes, there is mild evidence that the 
densities of these galaxies are lower. This is verified by considering the more carefully calculated mass at the half light radius based on the 
new Wolf et al. (2009) method, as described in § 6.4 and presented in Figure 13. 



Group dSphs extending down to the faintest systems 
(e.g., Mateo 1998; Grebel et al. 2003; Kirby et al. 2008), 
and the M31 satellites are found to track the brighter 
part of this relation remarkably well, e.g., from L ~ 10 5 
to 10 7 Lq. We illustrate this in the left panel of Fig- 
ure 11, where the darker encircled symbols represent the 
spectroscopically cleaned measurements from our present 
study. Such a strong correlation demonstrates that the 
gross chemical evolution histories of these satellites arc 
similar, and, if at all, similarly affected by dynamical 
interactions from their different host environments (see 
§6.5). As summarized in Table 3, the intrinsic metallic- 



ity spreads of the brightest satellites in the Milky Way 
and M31 are also similar to one another, and as large 
as 0.3 - 0.5 dex (the metallicity dispersion is indicated 
with horizontal error bars in Figure ll). 13 These simi- 
larities indicate that the stars in these galaxies may have 
resulted from an extended (early) star formation history, 
with the earlier generation of stars polluting the interstel- 
lar medium for subsequent star formation bursts. Yet, 
the trends of these overall processes must be very similar 

13 A horizontal error bar of 0.30 dex has been plotted for any 
dSphs that lack a metallicity dispersion measurement - See Table 3. 



18 



Kalirai et al. 



in the dSphs of the two hosts, and therefore minimally 
affected by external influences. As we note below, the 
gravitational potential of the dSphs in the two hosts also 
appears to be similar. Neither of these two findings, that 
there is a correlation between luminosity and metallic- 
ity or that the intrinsic abundance spreads are large, are 
true for globular clusters (Harris 1996). 

Interestingly, four of the least luminous M31 dSphs 
(And XI, XII, XIII, and XX) arc found to be signifi- 
cantly more metal-rich than expected given their lumi- 
nosities. These four systems, with —7.3 < My < —6.3, 
are found to be 0.5 - 1 dex more metal-rich than a rea- 
sonable extrapolation of both the M31 population and 
from the Milky Way dSphs with similar luminosity. One 
possible interpretation of this trend is that these satel- 
lites were initially much more massive than present, and 
have subsequently experienced tidal stripping. Such a 
scenario could possibly preserve the mean metallicity of 
the galaxies, while at the same time, reduce their lumi- 
nosities significantly. The presence of a positive radial 
abundance gradient (which is not known) towards the 
center of any of these dSphs would enhance this effect, 
as the more metal-poor stars in the outskirts would be 
preferentially stripped away (e.g., Chou et al. 2007). Un- 
fortunately, testing the tidal stripping possibility requires 
detailed kinematic data for many stars in these satellites, 
especially at larger radii. Such a data set does not yet 
exist in these low-luminosity systems. Similar to the M31 
sample, the faint Milky Way satellite Bootes II also devi- 
ates from the general trend, possibly suggesting a similar 
evolutionary history as outlined for these M31 satellites. 

6.3. Physical Sizes 

McConnachie & Irwin (2006) found the very interest- 
ing result that the sizes of M31's dSphs appear to be 
a factor of two larger than Milky Way dSphs of similar 
luminosity. The sample in their photometric study in- 
cluded some of M31's brightest satellites, And I, II, III, 
V, VI, and VII (e.g., down to My = -9.6). In Figure 11 
(middle-left panel), we extend this comparison by a fac- 
tor of 20 down to a luminosity of My — —6.3 by including 
all 19 M31 dSphs from the discovery papers summarized 
earlier. With the exception of the two newly found galax- 
ies And XIX (R e = 2.07 kpc - corrected from a geometric 
mean to an elliptical radius; McConnachie et al. 2008), 
and And XXI (R e = 875 pc; Martin et al. 2009), the 
bulk of the newly found dSphs in M31 have sizes com- 
parable to Milky Way counterparts of similar luminosity. 
We also note that the two outliers in the new data are all 
among the furthest galaxies along the line of sight, and 
have similar (or lower) luminosities than counterparts 
with smaller sizes. Therefore, the measurement of their 
sizes is more difficult given the fewer number of detected 
member stars. 

The overall size distribution of the two satellite sys- 
tems can be better seen on an expanded linear scale, 
as presented in Figure 12. The distributions are found 
to significantly overlap between L = 10 4 and 10 6 L@. 
In fact, at lower luminosities, it appears And XI, XII, 
XIII, and XX all have systematically smaller half light 
radii than similar luminosity Milky Way galaxies such as 
Bootes I, Hercules, Leo T, and Ursa Major I. We note 
that the uncertainty in R e in both of these figures is arbi- 
trarily taken to be 10% in cases where the error has not 
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Fig. 12.— The 2D elliptical half-light radii of Milky Way dSphs 
is compared to the M31 dSphs on a linear scale. In cases where 
the uncertainties in the half light radii were not directly measured, 
we have assumed it to be 10% in this plot. The general trend of 
dSphs in the two large hosts indicates that the most luminous M31 
satellites extend to larger sizes than the Milky Way dSphs (see 
also McConnachie & Irwin 2006). However, for the bulk of the 
lower luminosity systems (e.g., L < 10 6 I/©), the two distributions 
significantly overlap. 



been directly measured. From this overall comparison, 
we can say that less than one-third of the M31 satellites 
are consistently larger than Milky Way counterparts of 
the same luminosity, and that this result is mostly lim- 
ited to the brighter galaxies. This evidence therefore 
does not strongly support the claim that the Milky Way 
has a stronger tidal field that truncates its dSphs at a 
smaller radius relative to M31, unless such a dynamical 
effect could be limited to the brightest satellites. 

6.4. Dark Matter Masses 

The radial velocities from the SPLASH Survey pre- 
sented in this paper allow us to measure dynamical 
masses for roughly one-third of the entire known popula- 
tion of M31's dSphs. In Table 3, we have summarized the 
masses and mass-to- light ratios of these galaxies, based 
on the analysis presented in § 4 (i.e., using the Illingworth 
1976 method). We have also recalculated the masses of 
the Milky Way dSphs with resolved central velocity dis- 
persions using the same technique. 

As we noted earlier, these mass measurements follow 
from a simple scaling relation involving the measured 
sizes and velocity dispersions of these galaxies, and as- 
sume an isotropic sphere in which mass follows light. De- 
pending on the quality of the observations, several Milky 
Way dSphs allow modeling of the individual velocities 
to constrain the overall mass profile of the satellites (see 
e.g., Strigari et al. 2008b and Walker et al. 2009). Wolf 
et al. (2009) have presented a new prescription to cal- 
culate the masses of such systems out to the half light 
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radius, and show that these results agree well with with 
full dynamical modeling approach. We summarize the 
masses calculated from this technique, both for Milky 
Way and M31 satellites, in the last column of Table 3. As 
Wolf et al. (2009) show, the Illingworth (1976) approx- 
imation systematically underpredicts the mass of dSphs 
and so the new M1/2 are often similar, if not larger, than 
those values. This newer analysis is discussed in more 
detail below. 

In the two right panels of Figure 11, we illustrate the 
total measured masses and mass-to-light ratios of Milky 
Way dSphs and the six M31 dSphs that are the focus of 
this paper. Overall, the trend of increasing total mass 
of Milky Way dSphs with higher luminosity persists in 
the M31 dSphs. The slope of the relation between the 
two families of dSphs is also very similar. At the lu- 
minous end, our sample includes the two brightest M31 
dSphs known, And II and VII, both of which have a sim- 
ilar luminosity to the Galactic dSph Fornax. Formally, 
our measurements indicate that the total masses of these 
three dSphs is 6 - 9 x 10 7 M Q . At slightly lower lumi- 
nosities, we find that the M/Ly of And I (My = -11.8) 
is higher by a factor of two compared to the Milky Way 
dSph Leo I (My = —11.5), yet is lower by a factor of two 
compared to the luminous distant Local Group dSph Ce- 
tus (My = —11.3). For our three fainter satellites, we 
find that the M/Ly of And XIV is in nice agreement 
with the comparable Milky Way dSphs Draco, Canes Ve- 
natici I, and Leo T. However, for And III and And X, we 
find that the measured M/Ly is smaller by several fac- 
tors compared to most Milky Way dSphs of comparable 
brightness, such as Draco, Sextans, Ursa Minor, Canes 
Venatici I, and Leo T, and therefore these galaxies con- 
tain less dark matter than their Milky Way counterparts 
(see also Kalirai et al. 2009). 

For the luminous dSphs, the comparisons above 
demonstrate that, despite having a larger size than their 
Milky Way counterparts of the same luminosity, the 
masses of dSphs around each of the two hosts are sim- 
ilar. At lower luminosities, we see mild evidence that 
two of the six M31 satellites in our sample have a lower 
mass than their Milky Way counterparts. As these latter 
two satellites have half light radii that are similar, if not 
slightly larger, than Milky Way dSphs, this hints at an 
overall trend of the six M31 dSphs being less dense than 
Milky Way dSphs. 

Strigari et al. (2008b) have recently constructed mass 
profiles for the Milky Way dSph population and surpris- 
ingly found that, over four orders of magnitude in lumi- 
nosity, all of the satellites have approximately the same 
dark matter mass (~10 7 M ) within their central 300 pc. 
The interpretation of a common mass scale for Milky 
Way dSphs is still unclear, as discussed by Strigari et al. 
(2008b). For example, the common mass measurements 
may be indicative of a fundamental scale for dark mat- 
ter clustering. Correctly measuring whether such a scale 
extends to M31's dSphs, and if the absolute value is the 
same, will require both larger data sets (e.g., >50 accu- 
rate radial velocities per dSph) and more sophisticated 
modeling techniques. In the last column of Table 3, we 
have calculated the masses of each of the M31 dSphs and 
Milky Way satellites at the half light radius (M 1 / 2 ), using 
the new Wolf et al. (2009) technique (see §4.1). A com- 
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Fig. 13. — The masses (left) and mass-to-light ratios (right) of 
M31 (black) and Milky Way (grey) dSphs at the half light radius, 
calculated as described in Wolf et al. (2009). The comparison of 
dSphs in the two hosts indicates that the lower luminosity M31 
dSphs have dynamical masses and mass-to-light ratios that track 
the lower mass edge of the general relation for Milky Way satellites. 
Although the dynamical masses of these satellites are generally 
smaller than their Milky Way counterparts, they have half light 
radii that are similar or larger (see Figures 11 and 12). These 
galaxies are therefore less dense than the Milky Way dSphs. 



parison of M1/2 between the Milky Way and M31 dSphs, 
at the same luminosity, verifies the findings above. At 
the luminous end of the distribution, the M 1 / 2 are simi- 
lar despite the M31 dSphs being larger, while at lower lu- 
minosities, the M 1 / 2 are smaller for the M31 dSphs (e.g., 
And III, X, and XIV; see Figure 13). We will present 
the first results on the mass profiles of these M31 dSphs, 
and therefore the mass within a fixed physical radius, in 
Wolf et al. (2010, in prep.). 

6.5. The Bigger Picture: The Milky Way vs M31 

Although the Milky Way and M31 are both gener- 
ally characterized as large spiral galaxies, several differ- 
ences between the two systems have been pointed out 
in recent studies such as the SPLASH Survey. For ex- 
ample, surface brightness and abundance studies have 
shown that the metal-poor, stellar halo of M31 domi- 
nates the metal-rich inner spheroid at a radius of 20 - 
30 kpc, much farther out than the transition point in 
the Milky Way (Guhathakurta et al. 2005; Irwin et al. 
2005; Kalirai et al. 2006a; Chapman et al. 2006; Ibata 
et al. 2007; Koch et al. 2008). Additionally, the overall 
density of M31 at 20 kpc is much higher than the Milky 
Way (Reitzel, Guhathakurta, & Gould 1998), and M31 
contains a higher fraction of intermediate age stars than 
the Milky Way at these radii (Brown et al. 2007). Major 
differences among the population of satellites in the two 
systems have also long been known. The most luminous 
Milky Way satellites, the Small and Large Magellanic 
Clouds, are irregular type galaxies, of which M31 has 
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none. The most luminous M31 satellites are the dwarf el- 
lipticals NGC 205, NGC 187, and NGC 147, of which the 
Milky Way has none. The overall population of globular 
clusters in M31 outnumbers the Milky Way by several 
times. 

The differences between the Milky Way and M31 are 
likely tied to their past accretion histories. Deep, wide 
field imaging surveys of M31 demonstrate a large amount 
of substructure both in the inner and outer halo (Fergu- 
son et al. 2002; Ibata et al. 2007; McConnachie et al. 
2009). Simulations of galaxy formation in the ACDM 
paradigm suggest that the inner halos of massive galax- 
ies are built up quickly, and formed from just a few of 
the most massive, early accretion events (e.g., Diemand, 
Kuhlen, & Madau 2007). If the accretion history is re- 
sponsible for the present day differences between the two 
large galaxies (e.g., M31 having experienced more inter- 
actions - Hammer et al. 2007), then it may be the case 
that the surviving satellites suffered enhanced stripping. 
As we have discussed, the overall sizes of the M31 dSphs 
are systematically larger than their Milky Way counter- 
parts at the luminous end (see also McConnachie & Irwin 
2006), and the dynamical masses of M31's dSphs are ei- 
ther similar or lower than comparable luminosity Milky 
Way dSphs. However, in order to "puff" these systems up 
into less dense objects, the interactions would have had 
to strip the dark matter halos within the visible radii 
of the galaxies, in which case the stars should also be 
stripped. 

An alternate mechanism for losing stars through strip- 
ping has also just been presented by D'Onghia et al. 
(2009). They simulate the interactions of dwarf disk 
galaxies with larger systems and find that the encounters 
can excite gravitational resonances. They refer to this 
process as "resonant stripping" and show that the dwarf 
disk galaxies can transform to systems that look like 
present day dSphs (gas poor, dark matter dominated). 
Resonant stripping allows for up to 80% of the deeply 
embedded stellar component in satellite galaxies to be 
lost after a few Gyrs of interaction, while, at the same 
time, much less strongly affecting the surrounding dark 
matter halos. Therefore the mass-to-light ratios of these 
galaxies, and the sizes of their stellar distribution, can be 
altered depending on the details of the interaction his- 
tory (e.g., prograde vs retrograde orbits and the degree 
and timescale of the encounters). Unfortunately, such an 
effect is difficult, if not impossible, to measure from the 
limited angular spread of our radial velocities, which are 
not well extended in large numbers to the tidal radius 
of each satellite. A thorough investigation of the tidal 
interaction of these satellites with M31, and specifically 
relating such interactions to those of similar luminosity 
Milky Way dSphs, will also require an understanding of 
the orbits of these galaxies. Tangential velocities are cur- 
rently only available for a few of the most luminous Milky 
Way satellites, and for none of the M31 dSphs. 

To first order, we can investigate the projected distri- 
bution of M31's 19 dSphs relative to the center of the 
galaxy, as shown in Figure 14. The dSph distribution 
spans over 30 degrees in right ascension and declina- 
tion on this Figure (red circles). The concentric annuli 
are each separated by 30 kpc, and extend outward to 
180 kpc from M31. The small lines at the center of the 
figure mark the size of M31's optical disk. Also shown 
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Fig. 14. — The location of M31's 18 dSphs relative to the center 
of the galaxy is shown as filled red circles. The six dSphs com- 
prising this study are also indicated with larger open circles, and 
clearly span a wide range in distance from M31 and azimuth. As 
described in the text, the bias of satellites in the south results from 
incomleteness in the imaging surveys. 



are the locations of M33, and M31's 4 dwarf ellipticals 
(green triangles). First, this figure demonstrates that 
there is a clear absence of M31 satellites in the northern 
regions of the galaxy. The cause of this is an observa- 
tional bias related to the location of deep imaging surveys 
that have thus far only explored more southern latitudes 
(e.g., Ibata et al. 2007 and McConnachie et al. 2009). 
These surveys, e.g., PAandAS, are only now beginning to 
target other areas of the galaxy, and will surely uncover 
additional dSphs in M31, as will Pan-STARRS (Kaiser 
et al. 2002). The four recent discoveries of And XVIII, 
And XIX, And XX, and And XXI in the western hemi- 
sphere of M31 reflect the extention of the initial CFHT 
MegaCam survey (McConnachie et al. 2008; Martin et al. 
2009). 

The six satellites in our study are shown with enclosed 
open circles in Figure 14. These dSphs span an apprecia- 
ble range in radius and azimuth relative to M31. And I, 
III, and X sample intermediate radii in the halo, from R 
= 45 - 75 kpc, whereas And II, VII, and XIV are all lo- 
cated in the outskirts of M31's stellar halo, R > 150 kpc. 
In this context, one or more of the three inner satellites 
may represent close-in bodies that were captured long 
ago, whereas the outer group could be, in part, dSphs 
that are just now starting to fall into M31's potential. 
Evidence of such falling in galaxies is also supported by 
the radial velocity of And XIV as discussed earlier (see 
also Majewski et al. 2007), and also in the case of And XII 
(Chapman et al. 2007). However, both sets of galaxies, 
near and far from M31's center, include satellites that 
are systematically larger than Milky Way dSphs, as we 
showed in Figure 11. 

Interestingly, the three lower luminosity M31 satellites 
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And XI, XII, and XIII appear to be smaller than a Milky 
Way counterpart of similar luminosity (Bootes I), and are 
also located in a similar projected position in M31. Kine- 
matics of these galaxies will reveal if they are in fact mov- 
ing at similar radial velocities, possibly suggesting an in- 
falling group (e.g., see Li & Helmi 2008; D'Onghia 2008). 
The dark matter masses of these satellites, in relation to 
the general trend of Milky Way vs M31 dSphs, will help 
gauge the degree to which external dynamical influences 
may have played a role in shaping the internal kinemat- 
ics. More generally, the future of resolved spectroscopy 
of M31 dSphs is sure to shed light on both the intrinsic 
properties of the satellites and on the galaxy formation 
and evolution processes that shape these properties. 

7. CONCLUSIONS 

We present the first systematic large scale spectro- 
scopic survey of M31's dSphs as a part of the SPLASH 
Survey, a detailed study of And I, II, III, VII, X, and 
XIV. The fundamental properties of these galaxies are 
established from a sample of confirmed member stars, 
including the mean radial velocity, velocity dispersion, 
mean abundance, abundance spread, and total (stellar 
+ dark matter) mass. The kinematics indicate that 
the radial velocity and intrinsic velocity dispersion of 
the six satellites are v rad = —375.8 ± 1.4 km s _1 and 
o v = 10.6 ± 1.1 km s" 1 for And I, w rad = -193.6 ± 
1.0 km s" 1 and a v = 7.3 ± 0.8 km s" 1 for And II, 
Wrad = -345.6 ± 1.8 km s" 1 and a v = 4.7 ± 1.8 km s" 1 
for And III, w rad = -309.4 ± 2.3 km s" 1 and a v = 
9.7±1.6km s" 1 for And VII, v rad = -168. 3± 1.2 km s" 1 
and g v = 3.9 ± 1.2 km s _1 for And X, and v rad = 
-481.0 ± 1.2 km s" 1 and a v = 5.4 ± 1.3 km s" 1 for 
And XIV. 

Both photometric and spectroscopic metallicities are 
measured for all confirmed member stars in each galaxy, 
and the two independent measurements agree nicely with 
one another. We find that the six dSphs are metal-poor 
systems with ct m68.li iron &b undancB of [Fc/H] phot — 
-1.45 ± 0.04 (a [Fc/H] = 0.37) for And I, [Fe/H] phot = 
-1.64 ± 0.04 (<7[Fc/H] = 0.34) for And II, [Fe/H] phot = 
-1.78 ± 0.04 (o-[Fe/H] = 0.27) for And III, [Fe/H] phot = 
-1.40 ± 0.30 for And VII, [Fe/H] phot = -1.93 ± 0.11 
(o-[Fo/H] = 0.248) for And X, and [Fc/H] phot = -2.26 ± 
0.05 for And XIV (where two of the satellites lack a 
mctallicity dispersion measurement). 

The dynamical mass of each satellite is estimated to 
first order using the Illingworth (1976) formalism. More 
accurate mass measurements at the half-light radius of 
each galaxy are also calculated using the methods de- 
scribed in Wolf et al. (2009). As summarized in Table 2 
and 3, these yield half-light radius mass-to-light ratios of 
M 1/2 /L 1/2 = 31 ± 6 for And I, 13 ± 3 for And II, 19 
± 12 for And III, 7.7 ± 2.8 for And VII, 63 ± 40 for 
And X, and 102 ± 71 for And XIV. These calculations 
are found to agree nicely with masses measured through 
full modeling of the individual radial velocities assuming 
the dark matter follows a five parameter density profile. 

The comparison of the overall properties of M31's 
dSphs with Milky Way dSphs of similar luminosity in- 
dicates some interesting similarities and differences. For 
example, we find that the luminosity - mctallicity rela- 
tion is very similar between L ~ 10 5 - 10 7 L Q , suggest- 



ing that the chemical evolution histories of each group 
of dSphs is similar. The lowest luminosity M31 dSphs 
appear to deviate from the relation, possibly suggesting 
tidal stripping. The previous finding of McConnachie & 
Irwin (2006) indicating that the brightest M31 dSphs are 
systematically larger than Milky Way counterparts of the 
same luminosity does not persist cleanly to lower lumi- 
nosities. We find that the sizes of dSphs in the two hosts 
significantly overlap for L = 10 4 - 10 6 Lq, and that four 
of the faintest M31 dSphs are smaller than Milky Way 
satellites of similar luminosity. The first dynamical mass 
measurement comparisons between the M31 and Milky 
Way dSphs hints that the M31 satellites are systemati- 
cally less dense than Milky Way counterparts. 
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See references in §6.1 and Table 1 in Wolf et al. (2009). 
x Errors in r c (core radius), R e (2D elliptical half light radius), and r 1 / 2 (3D deprojected half light radius) are set to 10% (if not directly measured 
2 The total dynamical mass is calculated using the Illingworth (1976) method, as described in §4. This method systematically underpredicts the total 
3 The mass at the half light radius is calculated using Mj/ 2 = 3G _1 (o- 2 )r 1 / 2 (Wolf ct al. 2009), as described in §4.1. 
4 An error of 10% in distance was assumed. 

5 Given the uncertain measurement of the tip of the RGB, the distances to these satellites have been set to M31's mean distance, with a 10% uncertainty. 
6 An error of 0.5 mags was assumed for the integrated brightness of these galaxies. 
7 For this galaxy, we assume r c ~ Re, with a large 25% uncertainty. 



for the dynamical mass 
mass, as discussed in §§ 



estimates. 
4.1 and 6.4. 
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